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ABSTRACT 

We have analyzed archival Chandra and XMM-Newton data for two nearly complete homogeneously 
selected samples of type 2 Seyfert galaxies (Sy2s). These samples were selected based on intrinsic 
Active Galactic Nuclei (AGN) flux proxies: a mid-infrared (MIR) sample from the original IRAS 12/xm 
survey and an optical ([OIIIjA 5007 A flux limited) sample from the Sloan Digital Sky Survey (SDSS), 
providing a total of 45 Sy2s. As the MIR and [OIII] fluxes are largely unaffected by ACN obscuration, 
these samples can present an unbiased estimate of the Compton-thick (column density N/^ > 10^^ 
cm~^) subpopulation. We find that the majority of this combined sample is likely heavily obscured, 
as evidenced by the 2-10 keV X-ray attenuation (normalized by intrinsic flux diagnostics) and the 
large Fe Ka equivalent widths (several hundred eV to over 1 keV) . A wide range of these obscuration 
diagnostics is present, showing a continuum of column densities, rather than a clear segregation 
into Compton-thick and Compton-thin sub-populations. We find that in several instances, the fitted 
column densities severely undcr-represent the attenuation implied by these obscuration diagnostics, 
indicating that simple X-ray models may not always recover the intrinsic absorption. We compared 
AGN and host galaxy properties, such as intrinsic luminosity, central black hole mass, accretion 
rate, and star formation rate with obscuration diagnostics. No convincing evidence exists to link 
obscured sources with unique host galaxy populations from their less absorbed counterparts. Finally, 
we estimate that a majority of these Seyfert 2s will be detectable in the 10-40 keV range by the 
future NuSTAR mission, which would confirm whether these heavily absorbed sources are indeed 
Compton-thick. 
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1. INTRODUCTION 

A subset of galaxies are active, indicating that the cen- 
tral supermassive black hole (SMBH) is accreting ma- 
terial. Within such active galactic nuclei (AGN), this 
accretion disk is in turn surrounded by an obscuring 
medium of dust and gas, thought to have a toroidal ge- 
ometry (e.g. Antonucci 1993, Urry & Padovani 1995). In 
Type 1 AGN, the system is oriented such that the line of 
sight intercepts the opening of this "torus," exposing the 
accretion disk and broad load region (BLR) . Conversely, 
this central region is blocked in Type 2 AGN, as the sys- 
tem is inclined such that the line of sight is through the 
obscuring medium. In these obscured AGNs, the nar- 
row line region (NLR) is visible as well as scattered and 
reflected emission from the central engine. 

Previous studies have shown that obscured AGN con- 
stitute at least half of the local population (Risaliti 1999, 
Bassani 1999, Guainazzi et al. 2005). Obscuration can 
result from the putative torus or even the host galaxy 
where dust from nuclear star formation processes (e.g. 
Ballantyne 2008), extranuclear dust (Rigby et al. 2006) 
or dilution of AGN emission by the galaxy (Trump et al. 
2009) can attenuate optical signatures of AGN activity. 
However, in this study we focus on "toroidal" obscured 
sources (where the absorption is intrinsic to the AGN, on 
the scale of a few parsecs, enshrouding the BLR) since a 
substantial fraction of heavily obscured, Compton-thick 
(column density N// > 1.5x10^^ cm~^) AGN arc often 
invoked to explain the unresolved portions of the X-ray 



background at 30 keV (e.g. Gilli et al. 2007). An accu- 
rate census of these AGN is necessary to constrain X-ray 
background synthesis models (e.g. Treister et al. 2009) 
and studying their properties are crucial in understand- 
ing the full AGN population. Unfortunately, such ob- 
scured AGN are missed in 2-10 keV X-ray surveys as ab- 
sorption and Compton-scattering severely attenuate this 
X-ray emission. 

Selecting AGN samples based on intrinsic AGN flux 
proxies {P intrinsic), which are ideally unaffected by the 
amount of toroidal obscuration present, is therefore nec- 
essary to uncover Compton-thick AGN. Such diagnos- 
tics include emission lines which are primarily ionized 
by accretion disk photons and are formed in the NLR, 
making them not subject to torus obscuration, and in- 
clude the optical [OIII] 5007A line (e.g. Heckman et 
al. 2005) and the infrared [OIV] 25.89 /im line (e.g. 
Melendez et al. 2008, Diamond-Stanic et al. 2009, 
Rigby et al 2009). The obscuring medium absorbs con- 
tinuum photons and re-radiates them in the mid-infrared 
(MIR). This emission constitutes approximately 20% of 
the bolometric luminosity m . most type 1 and type 2 AGN 
(jSpinoglio fc Malkanlll989| ). making it another isotropic 
indicator of intrinsic AGN flux. Follow-up X-ray ob- 
servations can then reveal which sources arc potentially 
Compton-thick. 

Several studies have used infrared emission to locate 
such obscured AGN (Daddi et al. 2007, Fiore et al. 
2009). These studies have selected sources with infrared 
emission in excess of that attributable to star formation. 
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indicating the presence of an AGN. Using either Chan- 
dra and/or XMM-Newton spectra (Fiore et al. 2009) or 
stacked X-ray spectra for non-detections (Daddi et al. 
2007, Fiore et al. 2009), the column densities of these 
sources were estimated to constrain the Compton-thick 
fraction. However, these studies focus on high redshift 
sources {z = 0.7-2.5), rely on assumptions of the source 
spectrum shape (to convert count rate to flux for detec- 
tions or hardness ratio to flux for stacked spectra) and in 
the case of stacked sources, only characterize the aggre- 
gate population rather than individual sources. These 
analyses are useful in estimating a potential Compton- 
thick fraction at early times in the universe, but local 
sources generally have the advantage of X-ray spectra 
with higher signal-to-noise, where each source can be in- 
dividually analyzed and the spectra better characterized 
to more robustly constrain obscuration diagnostics. 

We have undertaken such an analysis on two homo- 
geneous samples of local (z < 0.15) Seyfert 2 galax- 
ies (Sy2s) selected based on intrinsic flux proxies: a 
MIR sample from the o riginal IRAS 12//m survey 
(jSpinoglio fc Malkanlll989D and an [OIII]-selected sample 
culled from the Sloan Digital Sky Survey ( LaMassa et ah] 
I2009D . The Chandra and/or XMM-Newton spectra of 
these sources were fit in a homogeneous and systematic 
manner. However, column densities derived from spec- 
tral fitting in the 2-10 keV band are highly model depen- 
dent and thus may not always refiect the intrinsic toroidal 
absorption. Other proxies are therefore necessary to 
identify potentially Compton-thick sources. As the 2-10 
keV X-ray emission is suppressed in absorbed sources, 
the ratio of this emission to intrinsic fiux indicators can 
probe the amount of obscuration. In Compton-thick 
sources, the ratio of the 2-10 kcV X-ray fiux (F2_iofcev) 
to ^intrinsic is about an order of magnitude or lower 
than what is observed in unobscured sources (e.g. Bas- 
sani et al. 1999, Cappi et al. 2006, Panessa et al. 
2006, Melendez et al. 2008). X-ray spectral signa- 
tures, most notably the equivalent width (EW) of the 
neutral Fe Ka line at 6.4 keV, can also aid in uncover- 
ing heavily obscured sources. As the EW is measured 
against a suppressed continuum, it rises with increasing 
column density, reaching values of several hundred eV 
to over 1 keV in Compton-thick sources (e.g. Levenson 
et al. 2002). Similar to LaMassa et al. 2009, we use 
^2-wkeV /^intrinsic and the Fc Ktt EW as obscuration 
diagnostics in this work. 

This paper is organized as follows: we describe the 
sample selection in Section 2 and the spectral fitting pro- 
cedures in Section 3. We then estimate the potential 
Compton-thick population using ancillary optical and 
IR data and discuss the various possible absorber ge- 
ometries revealed by our obscuration diagnostics. Utiliz- 
ing IR data to parametrize host galaxy characteristics, 
namely star formation processes, and AGN activity, such 
as intrinsic AGN luminosity, central black hole mass and 
Eddington ratio, we investigate whether Compton-thick 
sources trace a unique population. We also comment on 
the feasibility of detecting these sources in an upcoming 
hard X-ray (5-80 keV) mission, NuSTAR. 

2. SAMPLE SELECTION 

Our combined sample consists of 28 MIR selected and 
17 [OIII] selected Sy2s. The MIR sample is a subset 



of the 31 Sy2s from the original IRAS 12/^m survey 
(ISpinodio fc MalkanI 119891) Fl which was drawn from the 
IRAS Point Source Catalog (Version 2). The Sy2s in 
the 12/xm sample were selected via a weak color cut 
(i.e. 12/xm fiux less than the 60 or 100 /xm fiux) and 
is complete to a fiux-density limit of 0.3 Jy at 12/im, 
with latitude \b\ > 25° which av oids Galactic contamina- 
tion (jSpinoglio fc Malkanlll989() . Archival Chandra and 
XMM-Newton data exist for 28 of these 31 sources. 

The [OIII] -selected Sy2s were culled from the 
Main Galaxy sample (Strauss et al. 2002) in the 
Sloan Digital Sky Survey (S PSS) Data Release 4 
(lAdelman-McCarthv et al']l2006l) . Using the SDSS spec- 
tra, Sy2s were identified using the the diagnostic line 
ratio plot of [0III]/H/3 and [Nil] /Ha (BPT diagram, 
Baldwin et al. 1981) and the Kauffmann ct al. (2003) 
and Kewley et al. (2006) demarcations which distin- 
guish Sy2s, composite systems and star-forming galax- 
ies. The 20 local {z < 0.15) Sy2s with [OIII] 50071 fiux 
> 4 X 10"^'' erg s~^ cm"^ that lie within the AGN lo- 
cus of the BPT diagram were selected to comprise this 
sample. Of these 20 Sy2s, 2 had archival XMM-Newton 
observations and we were awarded XMM-Newton time 
for another 15 sources. The X-ray analysis for these 17 
[OIII] selected Sy2s was presented in LaMassa et al. 2009 
and is not replicated here; we utilize the results of that 
study (2-10 keV X-ray fiuxes, Fe Ka EWs, etc.) in this 
work. Though both original samples are complete, since 
X-ray data only exists for 28/31 and 17/20 sources from 
the 12/im and [OIII] samples respectively, our resultant 
sample for X-ray analysis is nearly complete. 

We note that this study has the advantage of analyz- 
ing samples of Sy2s selected via two different techniques 
which can mitigate biases from any individual selection 
criterion. For instance, dusty host galaxies can attenuate 
the optical emission lines used to identify AGN (and thus 
potentially miss galaxy obscured AGN), whereas MIR 
selection can isolate such sources. Star formation pro- 
cesses in the host galaxy can also enhance MIR, limit- 
ing its usefulness as an intrinsic indicator of pure AGN 
fiux. However, such effects are minimal in this study 
as all but two MIR identified Sy2s live in the AGN lo- 
cus of the BPT diagram and those two sources inhabit 
the composite (AGN and star-forming) locus and would 
not be optically identified as pure star-forming galaxies. 
We also note that many low-luminosity AGNs and some 
quasars lack the IR signature of a torus (Ho 2008, Hao 
et al. 2010, respectively), and thus MIR selection is not 
a useful tool for investigating such AGNs or their ob- 
scured counterparts. We explore the issues of selection 
effect biases further in Section 4.2. 

3. DATA ANALYSIS 

We analyzed the available archival Chandra and XMM- 
Newton observations fo r the 12/im sources with X Assist 
(|Ptak fc Griffiths! [200l . This program runs the appro- 
priate Science Analysis Systems (SAS) tasks to filter the 
data and clean for fiaring as well as extract spectra and 
associated response files for user-defined sources. Table[T] 
lists the X-ray observations used in this analysis, includ- 
ing the ObsIDs and net exposure times after filtering. 

^ Of the 32 Sy2s in the original sample, one was ro-classificd as 
a Syl (NOG 1097). 
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Twenty-five out of 28 sources were detected at tfic So- 
level or greater in the 0.5 - 10 keV band. One (NGC5193) 
was detected in the soft band (0.5 - 2 keV) and we were 
thus able to fit this part of the spectrum. We obtained 
upper 2-10 keV flux limits on this source and the two 
undetected sources (F08572-I-3915 and NGC 7590), dis- 
cussed in detail below. 

We used simple absorbed power law models to fit the 
spectra for the detected sources, which may not accu- 
rately represent the complex geometry of these systems. 
However, our main goal is to apply a systematic and ho- 
mogeneous ana l ysis o f the spectra in a similar manner as 
iLaMassa et al.l (|2OO90 to derive an observed X-ray fiux, 
and where possible, EW of the Fe Ka line. More exten- 
sive X-ray modeling of several sources have been inves- 
tigated in detail in the literature (e.g. see Brightman & 
Nandra 2010 for more detailed X-ray modeling of the ex- 
tended I2fj.ni sample) and we do not intend to replicate 
previously published work. In Appendix A, we discuss 
individual sources, compare our derived parameters with 
those quoted in the literature and comment on the im- 
pact more complex models have on such parameters. We 
find that in 18/23 sources, we recover consistent (within 
Icr) observed X-ray fluxes and Fc Ka EW values as more 
complex models. This work also represents the first anal- 
ysis for a handful of datasets (i.e. Chandra spectrum of 
IC 5063, XMM-Newton 2004 EPIC spectra of NGC 7172 
and XMM-Newton EPIC spectra of NGC 7674). 

3.1. Fitting Spectra from Multiple Observations 

Multiple observations for each source, as well as the 
spectra from the three XMM-Newton detectors (PN, 
MOSl and M0S2), were fit simultaneously with a con- 
stant multiplicative factor which was allowed to vary by 
~20% to account for calibration differences among de- 
tectors/observations. The remaining model parameters 
were initially tied together, with the residuals inspected 
to check for inconsistencies among observations. Differ- 
ences among XMM-Newton observations are interpreted 
as source variability, and were present in 4/28 sources 
(NGC 4388, NGC 5506, NGC 7172 and NGC 7582). 

Nine Sy2s had both Chandra and XMM-Newton 
archival data, with 8/9 having flux and/or spectral dis- 
crepancies between observations; only NGC 424 had con- 
sistent Chandra and XMM-Newton data. As Chandra 
has higher spatial resolution than XMM-Newton, it bet- 
ter isolates the central AGN. Differences in the spec- 
tra between the two observatories could thus be due to 
source variability, or extended emission from the host 
galaxy (e.g. X-ray binaries, thermal emission from hot 
gas, etc.) that XMM-Newton can not resolve from the 
AGN emission. To test if such differences were due to 
variability or contamination, wc extracted the Chandra 
source region to have the same size as the XMM-Newton 
region, ^20." If the best flt parameters and flux were 
consistent between the two datasets with the matched 
aperture extraction areas, we concluded that extended 
emission is likely contaminating the XMM-Newton ob- 
servation. If a discrepancy still existed, we interpreted 
this as source variability between observations. 

Five sources showed evidence of contamination from 
extended emission within the XMM-Newton aperture, 
i.e. matched aperture extraction between the Chan- 
dra and XMM-Newton observations resulted in consis- 



tent model parameters and fiux: NGC 1386, F05189- 
2524, NGC 3982, NGC 4501 and Mrk 463. For 3 of 
these sources (NGC 1386, F05189-2524 and Mrk 463), 
the best-fit parameters with the default Chandra ex- 
traction region were consistent with the XMM-Newton 
spectra, with the exception of the constant multiplica- 
tive factor which was lower in the Chandra observations 
(- 40 - 70% of XMM-Newton). We therefore fit the 
XMM-Newton and Chandra spectra simultaneously to 
constrain the Chandra parameters. However, we report 
the flux from the Chandra observations only in Table 
m as this isolates the central AGN. The spectra from 
the default Chandra extraction areas for the other two 
sources (NGC 3982 and NGC 4501) did not have consis- 
tent model parameters with the XMM-Newton spectra, 
likely due to X-ray binaries in the host galaxy affect- 
ing the spectral shape in the XMM-Newton data, so we 
therefore flt the Chandra spectra from the default ex- 
traction area independently and report these parameters 
in Table H 

Three Sy2s were variable between the two observa- 
tories: NGC 4388, Arp 220 and NGC 7582. Arp 220 
was flt simultaneously between the Chandra and XMM- 
Newton observations with only the absorption compo- 
nent fit independently for the Chandra spectrum. NGC 
4388 and NGC 7582 exhibited spectral variation between 
the Chandra and XMM-Newton observations and were 
therefore fit independently. Wc list the best fit parame- 
ters for the default Chandra extraction spectra and the 
XMM-Newton spectra separately in Table[2]for these two 
sources. 

3.2. Spectral Models 

We initially fit all spectra with an absorbed power law 
model. Most spectra (18/26) had an adequate number 
of detected photons to be grouped by a minimum of 15 
counts per bin without loss of spectral information and 
were thus analyzed with x^; the remaining 8 (NGC 3982, 
NGC 4501, TOLOLO 1238-364, NGC 4968, NGC 5135, 
NGC 5953, NGC 6890 and NGC 7130) were analyzed 
with the Cash statistic (C-stat, Cash 1979) and binned 
by 2-3 counts as XSpec handles slightly bi nned spec- 
tra b etter than unbinned when using C-stat (|Teng et al.l 
[2001 ). With the exception of 7 sources (NGC 1667, NGC 
3982, NGC 4501, TOLOLO 1238-364, NGC 4968, NGC 
5953 and Arp 220), a second power law component was 
needed to accommodate the data (i.e. phabsi *(powi -/- 
phabs2 *pow2)). The two power law indices (F) were tied 
together and the normalizations and absorption compo- 
nents were fit independently. Such a model represents 
a partial covering geometry with the first power law de- 
noting the soft scattered and/or reflected AGN contin- 
uum and the second component describing the absorbed 
transmitted emission. 

Residuals below 2 keV were present in many of the 
sources, suggesting emission in excess of the scattered 
AGN continuum. This excess is likely due to thermal 
emission from hot gas related to star f o rmati on processes, 
and consistent with ILaMassa et al.l (|2009l ). we used a 
thermal model (APEC in XSpec) with abundances flxed 
at solar to fit this emission. According to the f-test, ad- 
dition of this component improved the fit at greater than 
the 3(7 level over the best-fit single or double power law 
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model for 16/25 sourcesQ In Tabic [21 we present the 
best-fit parameters from the APEC plus power law mod- 
els, along with the values from the single absorbed 
power law fit, and where applicable, the double absorbed 
power law fit. We required a lower limit on the first ab- 
sorption component {Nh,i) to be equal to the Galactic 
absorption. In some cases, the best-fit absorption was 
equal to the Galactic Nh and we subsequently froze Nhs 
to the Galactic value for these sources. We were only able 
to obtain an upper limit on Nh i for three sources (Mrk 
463, NGC 6890 and NGC 7130),' as the lower error bound 
pegged at the Galactic absorption; the upper 90% limit 
is thus listed in Table [21 We also quote the 90% upper 
limit on kT for the six cases where the lower error on the 
temperature pegged at the limit of 0.1 keV (F05189-2524, 
NGC 3982, the Chandra observation of NGC 4388, NGC 
4968, NGC 5347 and the Chandra observation of NGC 
7582). We included Gaussian components to accommo- 
date the Fe Ka emission when present (see below) and 
additional Gaussian components for other emission fea- 
tures in NGC 1068 and NGC 7582 (see Appendix A). 
In Table [3 we list the best-fit parameters for the ab- 
sorbed single/double power law fit for NGC 5953 and 
the 9 sources which according to the f-test, are not sta- 
tistically significantly improved (> 3(t) by adding the 
APEC component and are therefore better described by 
the simpler single/double power law model (NGC 424, 
the Chandra observation of NGC 4388, NGC 4968, NGC 
5135, NGC 5347, NGC 6890, IC 5063, the Chandra ob- 
servation of NGC 7582 and NGC 7674). 

We list the observed 2-10 keV X-ray flux from these 
best-fit models in Table [H For the cases where addi- 
tion of the APEC model improved the fit, we excluded 
this component when deriving the X-ray flux. The flux 
was averaged among multiple observations when these 
observations were consistent. For variable sources, the 
flux is listed independently for each observation. For 
Arp 220, only the absorption varied between the XMM- 
Newton and Chandra observations, which had a negligi- 
ble impact on the flux. We therefore averaged the XMM- 
Newton and Chandra fluxes for this source. Wc note that 
NGC 7582 has a higher observed Chandra flux, compared 
to the XMM-Newton fluxes, despite the smaller Chan- 
dra spectral extraction area; aperture effects could con- 
tribute to the lower Chandra flux (compared with XMM- 
Newton) for NGC 4388. 

In Figure [H we plot the spectra with the best-flt mod- 
els. As many sources have multiple observations, we 
plot only one spectrum per observation, generally using 
the PN spectrum for XMM-Newton observations unless 
the MOS spectrum had better signal-to-noise. Though 
we report the flux of the Chandra observations only for 
NGC 1386, F05189-2524 and Mrk 463, we plot both the 
XMM-Newton and Chandra spectra to illustrate how the 
XMM-Newton spectra helped to constrain the fit. 

3.3. Pileup 

Bright X-ray sources can be susceptible to pileup which 
occurs when a CCD records two or more photons as a sin- 
gle event during the frame integration time. To test if 
this phenomenon affected our bright sources, we exam- 

^ Due to the marginal soft detection of NGC 5953, we did not 
fit this source with APEC. 



ined the pattern and observed distribution plots from the 
SAS task epatplot for XMM-Newton observations and the 
output of PIMM^ for Chandra observations. In a hand- 
ful of XMM-Newton observations (i.e. NGC 1068, NGC 
5506 and NGC 7172), one to two of the detectors exhib- 
ited evidence of pileup, but at least one of the detectors 
did not. The "piled" spectra were therefore disregarded 
from the fit without loss of information as we obtained 
one to two non-piled spectra per observation (see Ap- 
pendix A for details). As PIMMS uses simple models 
to test for the presence of pileup (e.g. single absorbed 
power laws whereas most of our sources needed a sec- 
ond power law component), we fit the Chandra spectra 
with evidence of pileup (i.e. IC 5063 and NGC 7582) in 
Sherpa, using the jdpilcup model and best-fit continuum 
model (with a Gaussian at the Fe Ka energy if neces- 
sary), to better constrain the pileup percentage. How- 
ever, we utilized the pileup model in XSpec (with a, the 
"grade migration" parameter, as the only free parameter) 
along with the best-fit models to derive the 2-10 keV 
fiux and Fe Ka EW, where the pileup component was re- 
moved before calculating these quantities. Wc note that 
the Sherpa and XSpec fits using their respective pileup 
models give consistent best-fit parameters and observed 
fluxes. 

3.4. Upper Limits 

Three sources were not detected within the 2-10 kcV 
range: F08572+3915, NGC 5953 and NGC 7590. NGC 
5953 was detected in the soft band (0.5 - 2 keV) and 
was therefore fit with an absorbed power law model. It 
was necessary to freeze the absorption to properly model 
the photon index. As the soft component generally re- 
sults from scattered/reflected AGN emission, the absorp- 
tion attenuating this component results from obscuration 
along the line of sight rather than intrinsic toroidal ab- 
sorption. In many cases in this study, such absorption is 
on the order of Galactic Nh or marginally higher, so we 
froze Nh to the Galactic value. From this flt in the soft 
band, we extrapolated an upper limit on the 2-10 kcV 
flux. 

F085724-3915 and NGC 7590 were not detected over 
the background in their ~15 ks Chandra and ~10 ks 
XMM-Newton observations, respectively. We therefore 
used a Bayesian approach to estimate an upper limit 
on the flux based on the total number of counts within 
the spectral extraction region and an assumed spectral 
shape for the AGN. We used a region size of ^2" for 
F08572-h3915 and ~7.5" for NGC 7590 (though XMM- 
Newton has lower resolution and the extraction region is 
generally ^^20", we constrained this region to a smaller 
size to exclude cont amination from a near by ultralumi- 
nous X-ray source (jColbert fc PtakI [2001 ). For NGC 
7590, wc coadded the MOS spectra together using the 
ftool addspec. We used the total detected and back- 
ground counts from these spectra to calculate a one-sided 
3(7 (i.e. 99.9% confidence level) upper limit on the num- 
ber of source counts. We then obtained an upper limit 
on the count rate by dividing this source count by the 
exposure time of the observation. Using an absorbed 
power law model, which included Galactic absorption, 
Compton-thick absorption (Nh = 1.5 x 10^"* cm~^, which 

^ http:/ /cxc. harvard. cdu/toolkit/pimms.jsp 
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is a conservative estimate as neither source was detected 
in X-rays) at the redshift of the source and a photon 
index of 1.8, we calculated the 2-10 keV flux that corre- 
sponds to the 3(7 upper limit on the count rate. These 
upper limits arc listed in Tabic |4l We note that applying 
this method to NGC 5953 results in a higher X-ray flux 
upper limit than extrapolating the spectral fit of the soft 
emission to higher energies, ~ 2 x 10~^^ erg s~^ cm~^ vs 
~ 5 X 10^^"* erg s^^ cm~^. We choose the latter value 
since this is based on the spectral information we have 
for this source. 

3.5. Fe Ka 

We used a Gaussian component to model the neutral 
Fe Ka emission. In many cases, this feature was evident 
when fitting the 0.5-8 keV spectrum and was included in 
the models mentioned above. For the sources where this 
line was not visible, we tested for its presence using the 
ZGAUSS model, freezing the energy at 6.4 keV and the 
width at 0.01 keV and inputting the galaxy's redshift. 
From this fit, we can derive either a detection or upper 
limit on the neutral Fe Ka fiux and possibly EW. For 
the sources that had both XMM-Newton and Chandra 
observations and had evidence of extended emission in 
the XMM-Newton field of view (i.e. NGC 1386, F05189- 
2524, NGG 3982, NGC 4501 and Mrk 463), we used only 
the Chandra spectrum to model the Fe Ka emission to 
isolate the AGN contribution. 

To better constrain the EW of the neutral Fe Ka line, 
we also fit the local continuum, from 3-4 keV to 8 keV, 
with a power law or double absorbed power law with 
an absorption component attenuating the second power 
law (when the spectral shape required this extra model). 
We then added a Gaussian or ZGAUSS component to 
this local continuum fit. The results of the global and 
local continuum fits to the neutral Fe Kaline are listed 
in Tabled In some cases (e.g. NGC 424, NGC 1386), 
the local fit better constrains the underlying continuum 
and therefore leads to a more reliable value for the EW. 
We use the EWs from the local fits in the subsequent 
analysis. 

Similar to LaMassa et al. 2009, we tested the signifi- 
cance of the Fe Ka EW detections by running simulations 
based on the power law(s) only component(s) of the local 
fit. We fit these simulated spectra with a Gaussian (or 
ZGAUSS) component to estimate the null hypothesis dis- 
tribution of line normalizations. Then the percentage of 
times that the observed line normalization exceeded the 
simulated line normalizations gives the statistical signif- 
icance of the line. 

4. DISCUSSION 

With the observed X-ray flux and Fe Ka EW con- 
strained, we can determine the distribution of the 
amount of 2-10 keV attenuation associated with the ob- 
scuring torus. As both the 12/im and [OIII] sample were 
selected on intrinsic AGN properties, such a percentage 
might represent an un biased estimate for the global AGN 
population. Similar to lLaMassa et all (|2009[ ). we also ex- 
plore if the fitted column densities agree with the proxies 
we use for AGN obscuration: if the emission is seen pri- 
marily via scattering and/or reflection, do the fitted N^ 
values recover the intrinsic absorption? The obscuration 
fiux diagnostics and Fe Ka EWs also provide clues as to 



the obscuration geometry in these sources. We compare 
host galaxy and AGN properties with Compton-thick di- 
agnostics to determine if sources with heavy absorption 
trace a unique populations from their less obscured coun- 
terparts. Finally, as higher energy (>10 keV) observa- 
tions are necessary to confirm a source as Compton-thick, 
we comment on the detectability of these Sy2s by NuS- 
TAR, an upcoming hard X-ray mission. 

4.1. Obscuration Diagnostics 

As fitted column densities are model dependent and 
could be unreliable, we use other proxies to investigate 
the amount of toroidal absorption in these systems, in- 
cluding the ratio of the observed X-ray flux to the in- 
herent AGN flux. We consider three diagnostics for in- 
trinsic AGN power (Fintrinsic)'- the [OIII] A 5007A line, 
the [OIV]25.89 ^m line and the mid-infrared (MIR) con- 
tinuum. The [OIII] and [OIV] lines are primarily ion- 
ized by the central engine, and as they form in the nar- 
row line region, are not subject to torus obscuration. 
The MIR emission results from the reprocessing of the 
AGN continuum by the dusty obscuring medium. We 
use the fiux at 13.5/im, averaged over a 3^m window, 
as Fmir since this region is free from strong emission 
li nes and absorpti o n feat ures. These fluxes are publ ished 
in lLaMassa et all (|2009[ ) and lLaMassa et all (|2010[) and 
are not replicated here. As these proxies are to first- 
order unaffected by the obscuring medium, whereas the 
2-10 keV X-ray flux is attenuated due to absorption and 
possibly Compton-scattering, the ratio of the X-ray flux 
to these tracers of intrinsic AGN power can probe the 
amount of obscuration present and has been used ex- 
tensively in previous studies (e.g. Bassani et al. 1999, 
Heckman et al. 2005, Cappi et al. 2006, Panessa et 
al. 2006, Melendez et al. 2008, LaMassa et al. 2009). 
We list the values of these obscuration diagnostic flux 
ratios in Table |6l There are, however, several limita- 
tions to using the [OIII] and MIR fluxes in tracing the 
intrinsic AGN flux: the [OIII] flux could be heavily af- 
fected by dust in the host galaxy and star formation pro- 
cesses can contaminate the MIR fiux (see lLaMassa et ahl 
(|20Tnh for a comparison between Fmi r/Fjoiii] between 
the two samples). In iLaMassa et all (|2010[ ). we noted 
that applying the standard R=3.1 reddening correction 
utilizing the Balmer decrement introduced errors that 
did not better recover the intrinsic [OIII] emission for the 
12/im sample, likely due to uncertainties in the II/3 mea- 
surements from the literature. Due to uncertainties in 
correcting the [OIII] and MIR fluxes for contamination, 
we use the observed parameters, with the caveat that 
these may not accurately probe intrinsic AGN emission 
for some sources. We discuss the implications of such 
biases below. 

We plot the distributions of F2-iokeV /F[oiii],obsi 
F2-iofcev/F[o/v] and F2-iokeV /Fmi in Figure [2] where 
the red dashed histogram represents the [OIII]-sample, 
the dark blue histogram denotes the non-variable 12/im 
sources and the cyan histogram reflects the variable 
12/im sources, using the average X-ray flux among the 
multiple observations for each source. A wide range of 
values is evident in all three plots. We compared our 
values with Syl sources, with the average flux ratio and 
spread delineated by the grey shaded regions in Fig- 
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ure [21 The Syl comp arison sample are culled from: a) 
iHeckman et al.l ()2005D (heterogeneous [OIII] -bright sam- 
ple, log(<F2 -wkev l^\ou i].obs >) = l-59±0.49 dex), 
b) iDiamond-Stanic et alj (j2009i ) (drawn from the re- 
vised Shapley-Ames catalog, log (<F2_iofeey/F[07y] >) 
= 1.92±0.60 dex ) and c iGandhi et al.l (IpOOl) fwhere 
F jvff B is calculated at 12.3/im with VISIR lLagage et aT 



- >VI I I , 

( 2004 ) observations of Sys selected from iLutz et al 



20041) and those with existing or planned hard (14-195 
keV) X-ray observations, log ( <F2-iofcey /FM/i?^ > ) = 
-0.34±0.30 dex). We note that iGandhi etall (|2009f ) re- 
port absorption-corrected X-ray luminosity whereas the 
other Syl comparison samples utilize the observed lumi- 
nosity. This correction shifts the F2_iofcev/FA//_R Syl 
ratios to higher values, though such a correction could 
be expected to be minimal for type 1 AGN which are 
thought to be largely unobscured. Also, not correct- 
ing [OIII] flux for reddening and MIR flux for starburst 
contamination could possibly result in obscuration diag- 
nostic ratios that are larger or smaller respectively, and 
though this affects several individual galaxies with large 
amounts of dust and/or greater star formation activity, 
no such systematic trends for the sample as a whole are 
evident. Yaqoob and Murphy (2010) have demonstrated 
that the ratio of V2-iokeV 1^ mir is more sensitive to 
the X-ray spectral slope and covering factor of the pu- 
tative torus, rather than column density, indicating that 
a low ratio does not necessarily imply a Compton-thick 
source. However, we find all three obscuration diagnos- 
tics to agree: the majority of Sy2s have ratios an order of 
magnitude or lower than their Syl counterparts, which 
may indicate Compton-thick absorption. 

This trend is further illustrated by Figure |3l which 
plots the observed X-ray luminosity as a function of in- 
trinsic AGN luminosity proxies, with the best-fit rela- 
tionship for Syls overplotted. Here, the red triangles 
represent the [OIII]-sample, the blue diamonds denote 
the non- variable 12/im sources and the cyan diamonds 
reflect the variable 12/im sources, with the individual X- 
ray fluxes (see Table H]) plotted for each variable source 
and connected by a solid line. The relationship for the 
Syl sources were calculated by multiple l inear regression 
(i.e. R EGR ESS routine in IDL) fo r the IHeckman et al.l 
(|2005f ) and IDiamond-Stanic et al.l ()2009l ) samples: for 
the MIR relationship, we uti lized the best-fit parame- 
ters from IGandhi et al.l (|2009[ ) for their Syl subsample. 
The majority of Sy2s lie well below the relations for Syls, 
demonstrating that these type 2 AGN have weaker ob- 
served X-ray emission. 

As the X-ray and optical and IR observations were not 
carried out simultaneously, it is possible that variability 
in the source could be responsible for the disagreements 
between the X-ray flux and intrinsic flux proxies. Such 
a scenario can be realized if the X-ray observations are 
made after the central source has "shut-off" (postulated 
to explain the discrepancy between the Type 1 optical 
spectrum yet r eprocessing-domin atcd X-ray spectrum for 
NGC 4051, see lMatt etall (|2003[ ) and references therein), 
or the converse, where optical observations are made 
during a sedentary state and X-ray observations cat ch 
the source in active state (e.g. iGuainazzi et al.l (|2005| )). 
Though we can not rule out variability as contributing 
to the discrepancy between the X-ray luminosity and in- 



trinsic AGN luminosity proxies for any individual source, 
such an effect can not be responsible for the overall trend 
in this sample. Variability in Syl samples contributes to 
the dispersion in lj2-iokeV /^isotropic ratios, yet they ex- 
hibit systematically higher X-ray luminosity (normalized 
by intrinsic AGN power) than Sy2s (Figures [2] and [3]) . 
This is confirmed by two-sample tests where we employed 
survival analysis (ASURV Rev 1.2, Isobe and Feigelson 
1990; LaValley, Isobe and Feigelson 1992; Feigelson and 
Nelson 1985 for univariate problems) to account for up- 
per limits in the X-ray flux. The Syl and Sy2 obscura- 
tion diagnostic ratios differ at a statistically significant 
level (< 1 X 10""' probability that they are drawn from the 
same parent population), which would not be expected 
if variability was the main driver for the discrepancy be- 
tween intrinsic AGN fiux and observed X-ray fiux. 

We have demonstrated that the majority of Sy2s in our 
samples arc under-luminous in X-ray emission as com- 
pared to Syls, but is this trend due to obscuration or 
inherent X-ray weakness? The EW of the neutral Fe 
Ka line can differentiate between these two possibilities 
and is thus another obscuration diagnostic. In heav- 
ily obscured sources, the AGN continuum is suppressed, 
whereas the Fe Ka line is viewed in refiection, leading to 
a large Fc Ka EW (several hundred cV to several keV, 
e.g. Ghisellini et al. 1994, Levenson et al. 2002). In 
Figure SI Fc Ka EW is plotted as a function of obscura- 
tion diagnostic ratios. A clear anti-correlation is present 
which is statistically significant according to survival 
analysis (Isobe et al. 1986 for bivariate problems): we ob- 
tain Spearman's p values of -0.648, -0.657 and -0.645 for 
Fe Ka EW vs. F2_iofcev/F[o///],ohs: F2-iofcey/F[ci/v] 
and V2-wkeV MIR, respectively. These best-fit cor- 
relations are overplotted for each relation in Figure ID 
The decrease of observed X-ray flux, normalized by in- 
trinsic AGN flux, with increasing Fe Ka EW indicates 
that obscuration is responsible for attenuating X-ray 
emission in these Sy2s. These results are consistent 
with the three-dimensional diagnostic diagram of Bas- 
sani et al. 1999 which shows a correlation between Fe 
Ka EW and column density which anti-correlates with 
F2-iofcev/F[o///],corr (whcrc F[o///]xorr IS the rcdenniug 
corrected [OIII] flux). 

Not only do a majority of this combined Sy2 sample 
exhibit trademarks of Compton-thick obscuration (an or- 



der of magnitude lower F2_iofcev/F 



isotropic 



ratios than 



Syls and large Fe Ka EW values), but a wide range 
of these diagnostic values are evident. No clear sepa- 
ration exists between Compton-thick and Compton-thin 
sub-populations. Also, though the diagnostic fiux ratios 
generally point to the same sources as having Compton- 
thick obscuration, not all three ratios agree for a hand- 
ful of sources (e.g. F05189-2524, NGC 5347, Arp 220, 
NGC 4388 and NGC 7582): some ratios indicate a 
Compton-thin source whereas others suggest Compton- 
thick. As the various intrinsic AGN indicators exhibit 
some scatter in inter-comparisons (see e.g. LaMassa et 



al. 2010), a spread in F2_iofcev/F. 



isotropic 



values is ex- 



pected. For F05189-2524, NGC 5347 and Arp 220, this 
discrepancy could be due to dust in the host galaxy af- 
fecting the [OIII] line, as mentioned above and/or large 
amounts of dust in the host galaxy boosting the MIR 
fiux. The 2005 XMM-Newton observation of NGC 7582 
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has a F2_iofcev/F[o///],o6s value marginally higher than 
the nominal Compton-thick/Compton-thin boundary, so 
the three flux ratio diagnostics may be considered to 
agree. However, the biases discussed previously in the 
observed [OIII] flux and MIR flux can not account for 
the disagreement of the diagnostic flux ratios in the 
Chandra and July XMM-Newton observations of NGC 
4388 and the 2001 XMM-Newton observation of NGC 
7582, where F2_iofcev/F[07y] point to the sources being 
Compton-thick at these stages, but the other ratios sug- 
gest a Compton-thin nature. Similarly, an Fe Ka EW 
of 1 keV is often cited as the nominal boundary for a 
Compton-thick source based on observations (e.g. Bas- 
sani et al. 1999, Comastri 2004, Levenson et al. 2006), 
yet NGC 1068, the archetype for a Compton-thick Sy2 
(|Matt et al.llT997l ). has a measured EW of 0.60;o;^^ keV 
(in agreement with Pounds & Vaughan 2006 but not 
Matt et al. 2004, see Appendix). Hence, though the di- 
agnostics presented here can help in uncovering the pos- 
sible Compton-thick nature of a type 2 AGN, nominal 
boundaries should be considered approximate, especially 
since a continuum of both diagnostic flux ratios and Fe 
Ka EWs are present. 



4.2. Implications for the Local AGN Population 

As both sub-samples were selected based on intrinsic 
AGN proxies and the majority is likely Compton-thick, 
this implies that heavily obscured sources could consti- 
tute most of the local AGN population. X-ray surveys 
in the 2-10 keV range, biased against these Compton- 
thick type 2 AGNs, would thus miss a significant por- 
tion of the population. Indeed, Heckman et al. 2005 
find that the luminosity function (which parametrizes 
the number of sources per luminosity per volume) for X- 
ray selected AGN is lower than the luminosity function 
for optically ([QUI]) selected sources. However, recent 
work (Trouille & Barger 2010, Georgantopoulos & Aky- 
las 2010) leads to the opposite conclusion, namely agree- 
ment between [QUI] and X-ray luminosity functions. As 
Georgantopoulos & Akylas (2010) point out, though the 
luminosity functions are similar, the selection techniques 
tend to find different objects, with [OIII] -selection fa- 
voring the X-ray weak sources. Hence, the number of 
sources per volume per luminosity may be comparable, 
but any one selection technique does not sample the full 
population. For instance, Yan et al. (2010) found that 
only 22% of their 288 optically selected AGNs are de- 
tected in the 200 ks Chandra Extended Groth Strip sur- 
vey, and they attribute the non-detection of the majority 
of the remaining sources to heavy toroidal obscuration. 
Conversely, X-ray selection can identify AGN that are 
categorized as star-forming galaxies by optical emission 
line diagnostics. For instance, Yan et al. (2010) note 
that about 20% of the X-ray sources identified as star- 
forming galaxies from optical emission lines have X-ray 
emission in excess of that explicable by star-formation, 
indicating the presence of an AGN. This finding is simi- 
lar to the results of Trouille & Barger 2010 who find that 
at least 20% of X-ray selected AGN in their sample are 
identified as star-forming according to optical diagnos- 
tics. Perhaps such competing biases work in concert to 
produce similar [OIII] and X-ray luminosity functions. 



4.3. Investigating Obscuration Ceometry 
4.3.1. Fitted Column Densities 

Here, we explore the relationship between obscuration 
diagnostics and the column densities (N//) derived from 
spectral fitting. In Figure \E\ we plot the fitted column 
densities as a function of F2-iokeV /^isotropic and Fe Ka 
EW for the 12/im sample (see LaMassa et al. 2009 for 
a discussion of fitted Nh for the [OIII] sample); as the 
[OIV] line is the least affected by the host galaxy contam- 
inations mentioned above, we use F^qiv] as F isotropic- 
With the exception of several sources, the fitted N// 
values approximately trace the degree of absorption im- 
plied by the obscuration diagnostics. However, a hand- 
ful of sources lay several orders of magnitude below this 
trend, and are labeled in Figure [5j This result is con- 
sistent with the findings of Cappi et al. (2006), where 
several Sy2s have fitted N^f values an order of magni- 
tude below that suggested by F2_iofcey /F[oiii].obs- Both 
F2-iokeV /Fisotropic and the Fe Ka EW diagnostics point 
to the same sources as being anomalous, with NGC 3982 
and NGC 4501 missing from the Fe Ka plot due to hav- 
ing an unconstrained EW or upper limit on the EW, 
respectively. All five of these sources required only a sin- 
gle power law model (with a thermal component in many 
cases) to adequately fit the spectrum. The low observed 
X-ray fluxes and high Fe Ka EW values indicate that the 
emission is primarily seen in scattering and/or reflection, 
rather than transmission through the obscuring medium. 
Hence such fitted values are likely associated with the 
line of sight absorption to the scattered/reflected com- 
ponent, suggesting that simple models of a foreground 
screen extincting the central source do not always re- 
cover the intrinsic absorption. 

Partial covering models, parametrized in this work by 
a double absorbed power law with the two photon indices 
tied together, can also misrepresent the inherent column 
density. For example, such a model fairly fit the spectra 
for NGC 1068 (x^=450.4 with 247 degrees of freedom), 
yet the best fit Nh was ~9x 10^^ cm~^ whereas the lower 
limit on this colum n density from hig her energy observa- 
tions is 10^^ cm-2 (jMatt et al.lll997| ). Though a partial 
covering model could more realistically represent the ge- 
ometry of the system, assuming a certain percentage of 
transmitted light through the obscuring medium with the 
rest scattered into the line of sight, it could be subject to 
the same limitations discussed above for single absorbed 
power law models. 

Published Njy distributions could potentially be bi- 
ased, skewed to lower values, though checks based on ob- 
scuration diagnostics can help mitigate this problem. For 
example, Akylas et al. (2006) analyzed the X-ray spec- 
tra for 359 sources from XMM-Newton and the Chandra 
Deep Field - South (CDFS), deriving intrinsic column 
densities from fitted Nh values though adopting a col- 
umn density of 5x lO^** cm~^ for the cases where F < 
1, a signature of Compton-thick obscuration. However, 
as Cappi et al. (2006) note, this criterion could indicate 
a Compton-thick source while the Fe Ka EW and fiux 
diagnostics suggest Compton-thin (e.g. NGC 4138 and 
NGC 4258) or vice versa (e.g. NGC 3079). Tozzi et al. 
(2006) use a reflection model (PEXRAV in XSpec) for 
Compton-thick sources in the CDFS, which are defined 
as those AGN with a better fit statistic using PEXRAV 
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than an absorbed power law model. However, as Mur- 
phy & Yaqoob point out (2009), such a model describes 
reflection off of an accretion disk, which is not physically 
relevant for the putative torus obscuration. Derived Nh 
values could then potentially be suspect for some sources. 
Other diagnostics are therefore crucial in checking the re- 
liability of fitted N// values. For example, Krumpe et al. 
(2008) find the ratio of X-ray to optical flux, as well as 
the non-detection of an Fe Ka line in the stacked spec- 
trum of 14 type II QSOs (AGN with intrinsic h^-wkeV > 
lO^''^ erg/s), to verify their distribution of moderately ab- 
sorbed, but not Compton-thick, column densities. 

4.3.2. Variable Sources 

It is intriguing to note that all X-ray variable sources 
in this study are on the high end of the obscura- 
tion flux diagnostics (see Figures [2] and [3]). These 
high 7 2-v:ikeV /^isotropic Aux ratios may indicate that 
the X-ray emission from the central source is seen 
directly. However, the high Fe Ka EW for the 
Chandra and July 2002 observations of NGC 4388 
{S).2'dtoll keV and 0.62t°;}° keV, respectively) and for 
the XMM-Newton observations of NGC 7582 (0.58lJ]:[!4 

keV and O.Slto'os keV) are higher than predicted for 
transmission-dominated spectra, where the EW with re- 
spect tojthe p rimary transmitted em ission is <0.18 keV 
(lMattl[200l . iPiconcelh et"aI1 (l200l propose a double 
absorption geometry to account for the variability in 
NGC 7582: a "thick" absorber which attenuates just the 
central source, attributed to the putative torus, and a 
"thin" absorber which enshrouds the primary and re- 
flected emission and is located externally to the torus. 
They postulate that this inner, "thick" absorber is inho- 
mogeneous, accounting for the observed X-ray variability. 
A similar scenario may be present for NGC 4388 and be 
responsible for both sources switching from transmitted- 
dominated to reflection-dominated states (or vice versa). 
The Fe Ka EWs for the two other variable sources, NGC 
5506 and NGC 7172, as well as the flux ratio diagnostics 
are consistent with Compton-thin sources, implying the 
central source is consistently viewed directly. 

4.4. Are Compton- Thick Sources Unique? 

Here we investigate whether Compton-thick sources 
differ from their Compton-thin counterparts in terms of 
host galaxy and AGN properties. In particular, we ex- 
amined whether systematic differences exist in intrinsic 
AGN power, Eddington ratio {Ijboiometncl^Eddington); 
central black hole mass {}Abh); the AGN contribution 
to the ionization held, and star formation activity. The 
results are summarized in Figures |6] through [11] and in 
Table [71 where we utilized survival analysis to calculate 
Spearman p values and the associated probabilities that 
the obscuration diagnostics are uncorrelated with host 
galaxy properties: P<0.05 indicates that the quantities 
are signiflcantly correlated (>2cr level). The values of the 
relevant hos t galaxy parame t ers us ed in this analysis are 
presented in iLaMassa et all (j201GD . 

To test whether Compton-thick sources have unique 
AGN properties, we searched for correlations between 
toroidal obscuration and intrinsic AGN luminosity, ac- 



cretion rate and central black hole mass (Mb_h0). As 
discussed previously, the [OIV] 25.89^m fine serves as 
a robust proxy of intrinsic AGN flux as it is mainly 
ionized by the central engine and not affected by host 
galaxy reddening as the [QUI] line is (e.g. Mclendez 
et al. 2008, Diamond-Stanic et al. 2009, Rigby et al 
2009). We therefore utihze L[o/y] as l^isotopic in Fig- 
ures [6| and [7| and Table [T] According to survival anal- 
ysis, a marginal statistically significant correlation ex- 
ists between Ijisotropic and two of the Compton-thick flux 
ratios [V2-wkev/'P[oiii],obs and F2_iofcey/F[o/y]), with 
a marginal significant anticorrelation between Lipotropic 
and Fe Ka EW. Figure [51 however, shows these depen- 
dencies to be weak with a wide scatter, especially consid- 
ering the error bars which can not be accommodated in 
the survival analysis calculation. We find no correlations 
between implied column density and accretion rate (us- 
ing Ljo/yj/Ms/f as a proxy for the Eddington ratio) and 
^BH (Figures [71 and [8]); survival analysis does indicate a 
weak significant relationship between Eddington param- 
eter and F2_iofeey/FA//i?,, but this is likely driven by the 
dependence on htQjyy As the weak correlation between 
luminosity and obscuration is tenuous at best, we con- 
clude that Compton-thick sources do not have systemat- 
ically different AGN properties from their less obscured 
counterparts. 

Could there be a relation between the obscuration 
shrouding the central engine and the large amounts of 
dust and gas necessary for starburst activity? As Lev- 
enson et al. (2004, 2005) point out, NGC 5135 and 
NGC 7130 (both members of the 12/im sample) are star- 
burst galaxies that likely harbor Compton-thick AGN. 
The combined [QUI] and 12/im samples provide us an 
opportunity to test if such a relation is generic. We use 
infrared quantities to illuminate the relative importance 
of starburst versus AGN activity: F[o/y]/F[jve//]7 which 
probes the hardness of the ionization field as F^qjv] 
is largely ionized by the AGN whereas [Nell] 12.81/im 
is excited by star formation processes; the EW of the 
17 ^tm polycyclic a, romatic hydrocarbon (PAH) feature 
(jCenzel et al.l Il998f) . which includes the emission fea- 
tures between 1 6.4-17.9 /xm ; and the MIR spectral in- 
dex a2o-30t,rM (|Deo et all [2009)1^ A higher value of 
F[o/v]/F[iVe//] indicates the dominance of AGN activity 
whereas larger PAH EW at 17/im and 0120-30^1171 values 
denote higher levels of starburst activity. As Figures [51 - 
[TT1 and Table [71 illustrate, a correlation between column 
density and hardness of incident ionization field/star for- 
mation activity do not exist. These results suggest that 
the gas responsible for starburst processes likely origi- 
nates in regions of the galaxy not associated with the 
putative torus, and similarly that gas from the interstel- 
lar medium does not contribute significantly to toroidal 
AGN obscuration in hard (2-10 keV) X-rays. 

Mflg measured us ing velocity dispersion and the M-ct relation 
IITremaiiie et al. I [20021) . See LaMassa et al. 2010 for literature 
references to Mg^f for the 12/jm sample; velocity dispersions for 
the [OIII]-sample were derived from SDSS. 

^ "Ai-Aa = logifxj f\2)/loa{^i/W 

^ We note that we only have these data for 27 of the 28 12/^m 
sources presented in this work as NGC 1068 had saturated low- 
resolution Spitzer data. We were therefore unable to obtain a PAH 
17 /im EW value or a20-30(im value for this source. 
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4.5. NuSTAR: Detection at Higher Energies 

In order to confirm a source as Compton-thiek, ob- 
servations at higher energies (>10 keV) are necessary. 
The spectral characteristic of heavily obscured AGN is 
the so-called "Compton-hump" , a peak in the spectrum 
between 20-30 keV which is caused by the competing 
effects of absorption on the low-energy end and Comp- 
ton down-scattering on the high energy range. The Nu- 
clear Spectroscopic Telescope Array (NuSTAR), to be 
launched in 2012, is sensitive in the 5-80 keV band, and 
could thus confirm our obscured candidates as Compton- 
thick sources, if they are detected. 

To test if these sources would be detectable by 
NuSTAR, we simulated higher energy spectra, 
using the XSpec command fakeit, based on the 
best-fit model and the associated response and 
backgroun d files provided by the NuSTAR team 



(http://www.nustar.caltcch.edu/for-astronomers/simulationsl. 
For the three non-detections in the 12/im sample, we 
simulated spectra using a model that takes into account 
Compton scattering assuming a spheroidal obscuration 
geometry (PLCABS in XSpec), with Nh=10'^'^ cm'^, 
r=1.8 and the maximum number of scatterings set to 
5; the normalization was adjusted such that model 2-10 
keV flux equaled the upper limits we calculated via 
Bayesian analysis. Using the simulated observed source 
and background count rate, we estimated the exposure 
time necessary for a source to be detected at the 5a level 
over the background. We find that all but five sources 
from the Ufirn sample (NGC 1386, NGC 1667, Tololo 
1238-364, NGC 4968 and NGC 6890) and four from 
the [Olllj-sample (2MASX 08035923-1-2345201, 2MASX 
J10181928-I-3722419, 2MASX J13463217-I-6423247 and 
NGC 5695) will be detected with exposure times less 
than 100 ks (see Appendix). However, though our 
simulations indicate the three non-detected 2-10 keV 
sources will be observable by NuSTAR, this is an 
opti mistic estimate, and sh ould be treated with caution. 

In iLaMassa et all ()2010[ ). we noted that the majority 
of these Sy2s are undetected by the Swift BAT Surveys, 
indicating that these sources are heavily absorbed. How- 
ever, as NuSTAR probes to a much deeper flux level 
in a million second observation than the BAT surveys 
(~2xl0~" erg s^^ cm~^ vs. the limiting BAT flux of 
~ 3.1 X 10""'^^ erg s~^ cm~^), the majority of these heav- 
ily obscured sources will likely be detected if observed by 
NuSTAR. 

5. CONCLUSIONS 

We have analyzed archival Chandra and XMM-Newton 
observations for two nearly complete homogeneous sam- 
ples of Sy2 galaxies: one selected from the SDSS on the 
basis of observed [OIII] flux and a MIR sample from 
the original IRAS 12/xm sample. The combined sample 
provided us with 45 Sy2s with existing Chandra and/or 
XMM-Newton data. Of these, three were not detected 
above the background (F08572+3915, NGC 5953 and 
NGC 7590) and four exhibited evidence of variability 
among multiple X-ray observations (NGC 4388, NGC 
5506, NGC 7172 and NGC 7582). 

We probed the amount of absorption present in these 
sources by comparing the 2-10 keV X-ray flux with opti- 
cal and MIR proxies of intrinsic AGN power {F intrinsic- 



the fluxes of the [OIII]A 5007 A and [OIV]25.89 /^m emis- 
sion lines and the MIR continuum flux at 13.5 /xm) and 
by investigating X-ray spectral signatures of obscuration 
(i.e. Fe Ka EW). We compared such obscuration diag- 
nostics with fltted column densities and explored the im- 
plications of these diagnostics on the AGN geometry. We 
also investigated whether a connection exists between the 
column density of the obscuring medium and host galaxy 
characteristics. Our results are summarized as follows: 



1. The majority of the combined sample has 
F2-iokeV /^intrinsic valucs an ordcr of magnitude 
or lower than the mean values for Syls. The 
statistically significant anti-correlation between 
F2-iokeV /^intrinsic and Fe Ka EW indicates that 
these lower diagnostic flux ratios are due to ob- 
scuration rather than inherent X-ray weakness in 
Sy2s. Thus a majority of these sources are poten- 
tially Compton-thick, consistent with the results of 
previous studies (e.g. Risaliti 1999). 

2. A wide range of obscuration diagnostic values are 
present, indicating a continuum of column densi- 
ties and/or inclination angles, rather than a clear 
segregation into Compton-thick and Compton-thin 
sub-populations. Though the diagnostics do gen- 
erally point to the same sources as likely heavily 
absorbed, disagreement does exist for a handful of 
Sy2s. Such a discrepancy is to be expected based 
on the various biases affecting the observed intrin- 
sic flux proxies and the inherent spread in such 
isotropic flux indicators (e.g. see ILaMassa et al.l 
(I2010D ). Hence, nominal Compton-thick bound- 
aries should be considered approximate. 

3. Thoug h recent work (iGeorgantopoulos fc AkvlasI 
(pOlOl ). iTrouille fc Barged (|2010l )) shows the lu- 
minosity functions for X-ray selected and [OIII]- 
selected AGN to be consistent, the various selection 
techniques favor differ classes of objects. Heav- 
ily obscured sources, present in optically selected 
samples, are missing from 2-10 keV X-ray sam- 
ples. Sample selection based on intrinsic flux prox- 
ies are therefore necessary to include the Compton- 
thick population, especially since highly absorbed 
sources constitute the majority of our homoge- 
neously selected samples. 

4. Though fltted column densities generally tend to 
trace the absorption implied by obscuration diag- 
nostics, several glaring inconsistencies are present. 
Such discrepancies are most extreme when the hard 
X-ray spectrum is best fit by a single absorbed 
power law. implying that the simple geometry of 
a foreground screen attenuating the central source 
does not recover the intrinsic absorption. This 
could result from scattering and/or reflected emis- 
sion dominating over the transmitted continuum, 
where the fltted column density reflects line of 
sight absorption rather than the obscuration en- 
shrouding the AGN. Such a result indicates that 
published Nh distributions derived from single ab- 
sorbed power law models can be similarly biased, 
systematically under-representing the intrinsic col- 
umn density of type 2 AGN. Other diagnostics are 
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therefore crucial in checking the vahdity of fitted 
column densities. 

5. The X-ray variable Sy2s populate the less ob- 
scured range of the flux ratio obscuration diagnos- 
tics. Two of these sources (NGC 4388 and NGC 
7582) do show evidence of switching to a reflection- 
dominated state, as indicated by the change in 
the Fe Ka EWs. As Piconcelh et al. (2007) 
suggest, this change could reflect an inhomoge- 
neous thick absorber covering the central source, 
with a thin absorber attenuating both the reflected 
and transmitted emission. The other two variable 
sources (NGC 5506 and NGG 7172) show signs of 
Compton-thin absorption, suggesting that the cen- 
tral source is viewed directly. 

6. We do not find compelling evidence that Compton- 
thick sources have unique AGN properties (intrin- 
sic AGN luminosity, accretion rate and central 
black hole mass) or star formation activity. Though 
three out of four obscuration diagnostics are signif- 
icantly correlated with intrinsic AGN luminosity, 
the significance is marginal and the relationships 
display a wide scatter. Evidence linking more ob- 
scured sources to more luminous central engines 
is therefore tenuous at best. No correlation exists 
between toroidal AGN obscuration and the relative 
amount of ionization due to the central engine com- 
pared to star formation processes {F[oiv]/^[Neii]-, 
EW of the 17^m PAH feature, 020-30^™) and AGN 
absorption. Though several starburst galaxies do 
seem to host Compton-thick AGN (e.g. Levenson 
et al. 2004, 2005), such a relation is not present 
globally. Hence, we conclude that the gas respon- 
sible for star formation processes is not associated 
with the toroidal obscuration hiding the central en- 
gine. 

7. Based on simulated high-energy (10-40 keV) spec- 
tra using the best-fit modeling of the 2-10 keV spec- 
tra, we estimate that the majority of this sample 
(36 out of 45) will be detected if observed by NuS- 
TAR. The more heavily obscured sources which 
have not been detected by BAT surveys could likely 
be identified by NuSTAR as this future mission will 
probe to lower flux levels (~2xl0~^^ erg s^^ cm^^ 
vs. ~ 3.1 X 10~^^ erg s~^ cm~^). These observa- 
tions would confirm whether the heavily absorbed 
sources are indeed Compton-thick. 
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Figure 1. X-ray spectra with best-fit models. 
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Figure 2. Histograms showing the distribution of obscuration diagnostic ratios {F2-l0keV /^iaotropic)- Dark blue histogram rep- 
resents the non- variable I2fj,m sources, cyan reflects the X-ray variable 12/im sources (X-ray fluxes are averaged for each so urce), 
red denotes the [OIII]-sample and the gray shad ed region illustrates the aver age value for Syls from a) IHeckman et al.l II2005I ). log 
(<F2-io*:ev/F[o///l,of,s >) =1.59±0.49 dex, b) IDiamond-Stanic et al.l l|2009l '). log (<F2_iofcev/F[o/v] >) =1.92±0.60 dex and c) 
IGandhi et al.l I I2009I) . log {<F2-iOkeV /^MIR >) =-0.34 ±0.30 dex. The left facing arrows represent X-ray upper limits in a), b) and 
c) and right facing arrows illustrate the [OIV] upper limits in b); these values are not included in the histogram. 
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Figure 3. X-ray luminosity versus proxies of intrinsic AGN luminosity. Blue diamonds represent the non-variable 12;im sources, cyan 
diamonds illustrate the variable 12^m sources, with the observed flux values from Table|4]for each source connected by a straight line, and 
red triangles denote the [OIII] -sources. Error bars for the variable sources represent the upper error on the maximum X-ray flux and lower 
error on the minimum X-ray flux. Variable sources NGC 5506 and NGC 7172 do not have error bars plotted as they are smaller than the 
sy mbol size. The dashed line rep resents the relationship for Syls from a) IHeckman et al. ( 200j), s l ope = 1.4 dex with intercept -0.15 dex 
, b) IDiamond-Stanic et aLl II2009I ). slope = 0.86 dex with intercept = 7.53 dex and c) [Gandhi et al.l I I2009I) . slope = 0.85 dex with intercept 
6.27 dex. In all cases, the majority of the Sy2s are below this relationship, illustrating that Sy2s have weaker X-ray emission than their 
Syl counterparts. 
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Figure 4. Fe Ka EW as a function of obscuration diagnostic ratios. Color coding is similar to \3\ though a dashed line is used for 
variable source NGC 7582 to avoid confusion with other variable sources having similar values. The statistically significant anti-correlations 
among all three relationships (Spearman's p= -0.647, -0.657, -0.645, respectively, calculated from survival analysis) indicate obscuration 
is primarily responsible for X-ray attenuation. In b), the two sources with lower limits on F2-I0ke.v /^[OIV] ^''s plotted for illustrative 
purposes and were not included in the survival analysis calculations. The fitted relationships are a) slope = -0.36 ±0.07 dex with cr=0.37 
dcx and intercept of -0.05 dex, b) slope = -0.47 ±0.08 dex with (7=0.34 dex and intercept of 0.002 dex and c) slope = -0.41 ± 0.06 dex 
with (7=0.32 dex and intercept of -1.00 dex. 
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Figure 5. Fitted N// as a function of obscuration diagnostics for the sample. The sources without error bars either had the best-fit 

absorption equal to the Galactic value, and therefore frozen at this value during fitting, or had J^h error bars smaller than the symbol size. 
The dashed lines indicate the boundary for a Compton-thick column density (N^ > 1.5 X 10^* cm"'^) and the dashed-dotted line indicates 
nominal Compton-thick boundaries based on obscuration diagnostics (log {F2—ioke.v /^[OIV]) 0-9 dex, an order of magnitude less than 
the average value for Syls, and Fe Ka EW > 1 keV). Sources that are likely heavily obscured according to the obscuration diagnostics, yet 
have low fitted column densities, are labeled. Color coding same as Figure |4l 
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Figure 6. Obscuration diagnostics vs. intrinsic AGN luminosity, parametrized by l^ioiv] Tiie lower limits on F2~wkeV /^[OIV] '^.re 
displayed for illustrative purposes and not included in the survival analysis calculation. Survival analysis indicates a marginal statistically 
significant correlation among three of these relationships (p=0.273, 0.185, 0.349 and -0.302, respectively), however a wide scatter is evident. 
Color coding same as Figure |3] 




Figure 7. Obscuration diagnostics vs. Eddington ratio. The lower limits on F2-iOkeV /^[OIV] ^.re displayed for illustrative purposes and 
not included in the survival analysis. With the exception of F2—l0keV MIR> which survival analysis suggests is marginally significantly 
correlated with Eddington parameter, no trends are present: p=0.110, 0.056, 0.280 and -0.219, respectively. Color coding same as Figure 

m 




Figure 8. Obscuration diagnostics vs. Ms//. The lower limits on F2~iOkeV /^[OIV] displayed for illustrative purposes and not 
included in the survival analysis. No statistically significant trends are present: p=0.148, 0.115, -0.012 and -0.234, respectively. Color 
coding same as Figure |3] 
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Figure 9. Obscuration diagnostics vs. F[oiv]/^[Nell]^ ^ proxy for the relative strength of the ionizing continuum from the AGN versus 
starburst activity. The lower limits on F[oiv]/F[Nell] s-^d F2_iofce\^/F[o/V'] ^^'^ displayed for illustrative purposes and not included in 
the survival analysis. No statistically significant trends are apparent: p=0.036, 0.032, 0.269 and -0.130, respectively. Color coding same as 
Figure [3] 
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Figure 10. Obscuration diagnostics vs. the EW of the PAH 17 fim feature, which parametrizes star formation rate. The lower hmits on 
F2— l0fcev/F[O7V'] S'fc displayed for illustrative purposes and not included in the survival analysis. No statistically significant trends are 
apparent: p=0.135, 0.062, -0.055 and -0.192, respectively. Color coding same as Figure [3l 
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Figure 11. Obscuration diagnostics vs. alpha2o-30Mm i which parametrizes star formation rate. The lower limits on F2-l0keV /^[OIV] 
are displayed for illustrative purposes and not included in the survival analysis. No statistically significant trends are apparent: p=C).272, 
0.014, 0.060 and -0.235. Color coding same as Figure [3] 
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Table 1 

Sample and Observation Log 



Galaxy 


Distance 
Mpc^ 


Observatory 


Observation Start Date 
UT 


ObsID 


— 

Exposure Time 
MOS1/MOS2/PN2 
ks 


JNCjO U424 


51.2 


AJVIM 


2001 Dec 12 


00029242o01 


7.6/7.0/5.0 






Chandra 


20U2 leb 4 


03146 


9.1 


JNCjO 1Uo8 


16.9 


AJVIM 


2000 Jul 29 


0111200101 


38.7/35.0/35.3 




16.9 


XMM 


2000 Jul 30 


0111200201 


37.8/35.0/32.2 


NGC 1144 


120.8 


AT" TV T TV T 

XMM 


2006 Jan 28 


0312190401 


11.6/11.6/10.0 


NGC 1320 


38.3 


■\7" TV T TV T 

XMM 


2006 Aug 6 


0405240201 


16.8/16.8/13.7 


NGC 1386 


12.7 


■\ 7" TV T TV 

XMM 


2002 Dec 29 


0140950201 


17.1/17.1/15.1 






Chandra 


2003 Nov 19 


04076 


19.6 


NGC 1667 


64.1 


XMM 


2004 Sep 20 


0200660401 


10.0/10.1/8.1 


F05189-2524 


187.7 


XMM 


2001 Mar 17 


0085640101 


10.7/10.6/7.6 






Chandra 


2001 Oct 30 


02034 


18.7 






Chandra 


2002 Jan 30 


03432 


14.9 


F08572+3915 


256.0 


Chandra 


2006 Jan 26 


06862 


14.9 


NGC 3982 


16.9 


XMM 


2004 Jun 15 


0204651201 


11.5/11.5/9.7 






Chandra 


2004 Jan 3 


04845 


9.2 


NGC 4388 


34.0 


Chandra 


2001 Jun 8 


01619 


20.0 






XMM 


2002 Dec 12 


0110930701 


11.7/11.7/7.8 






XMM 


2002 Jul 7 


0110930301 


9.0/9.2/2.8 


NGC 4501 


34.0 


XMM 


2001 Dec 4 


0112550801 


13.4/13.4/2.9 






Chandra 


2002 Dec 9 


02922 


17.9 


TOLOLO 1238-364 


46.9 


Chandra 


2004 Mar 7 


04844 


8.7 


NGC 4968 


42.6 


■\7" TV T A T 

XMM 


2001 Jan 5 


0002940101 


7.3/7.3/4.9 






■\ 7" TV T TV 

XMM 


2004 Jul 5 


0200660201 


4.5/4.7/5.2 


M-3-34-64 


72.7 


XMM 


2005 Jan 24 


0206580101 


44.6/44.6/42.9 


JNGC 5135 


59.8 


Chandra 


2001 Sep 4 


02187 


29.3 


ISILjC 5194 


8.5 


Chandra 


onnrv T on 

zuUu Jun zv) 


00354 


14.9 






Chandra 


2001 Jun 2S 


01622 


26.8 






Chandra 


zUOo Aug ( 


nonQO 


47. y 


JNCjC oo47 


34.0 


Chandra 


2004 Jun 5 


04867 


36.9 


iVlrk; 4do 


oi n /I 


AiVllVi 


zUUi Uec iz 


nnn/1 /i ni om 


zO.U/ZO.U/zo.4 






Chandra 


zUU4 Jun il 


04913 


49.3 


INLtO OOUd 


25.5 


AMM 


zUUl reb z 


nm '31 /I m m 
OUioi40iOi 


1 ( .8/1 i .8/14.3 






XMM 


2002 Jan 9 


nni oi /inom 

0013140201 


10 0/10 o/in 

13.2/13.2/10.6 






XMM 


2004 Jul 11 


nom oononi 
O2OI00O2OI 


01 /OI /OI 1 
21.0/21.0/21.1 






AMM 


2004 Jul 14 


nom OQnQm 
O2OI00O0OI 


on /on /in '7 
20.2/20.2/19. { 






"VTV TTV J 

XMM 


Of\r\ ATI o O 

2004 Jul 22 


rvorvi oorvnirvi 

0201830401 


in /in /in n 

19.0/19.0/19.9 






■\7" TV T TV T 

XMM 


2004 Aug 7 


0201830501 


20.2/20.2/20.0 






"VTV TTV J 

XMM 


2008 Jul 27 


0554170201 


or* o/oo n/nn -I 

85.2/88.0/90.4 






"VTV TTV J 

XMM 


2009 Jan 2 


K A ~i 101 

0554170101 


75.1/76.0/87.0 


iNLrU oyoo 


29.7 


Chandra 


zU02 Uec 12 


04023 


4.7 


Arp 220 


77.1 


■\ 7" TV T A J" 

XMM 


2002 Aug 11 


0101640801 


13.6/13.6/11.8 






■\7" TV T A T 

XMM 


2003 Jan 15 


0101640901 


14.6/14.6/9.3 






XMM 


2005 Jan 14 


0205510201 


O T / o o / r\ T 

8.7/8.2/0.7 






"VTV TA T 

XMM 


2005 reb 19 


nonr'r'i n,ini 

0205510401 


OI / a ^ / A o 

8.1/8.3/4.3 






Chandra 


2000 Jun D 


nn n 

00869 


56.5 


IN vj v_,' uoyu 


34.0 


XMM 


zi\juo ocp zy 


(T^m 1 "^1 nm 

UOU J- ± ij ±UU -L 




IC 5063 


46.9 


Chandra 


2007 Jun 15 


07878 


34.1 


NGC 7130 


68.4 


Chandra 


2001 Oct 23 


02188 


38.6 


NGC 7172 


38.3 


XMM 


2002 Nov 18 


0147920601 


13.6/13.6/12.0 






XMM 


2004 Nov 11 


0202860101 


50.8/50.9/36.0 






XMM 


2007 Apr 4 


0414580101 


48.9/48.8/31.7 


NGC 7582 


21.2 


Chandra 


2000 Oct 14 


00436 


10.5 






Chandra 


2000 Oct 15 


02319 


5.9 






XMM 


2001 May 25 


0112310201 


22.6/22.6/19.6 






XMM 


2005 Apr 29 


0204610101 


80.2/79.7/71.8 


NGC 7590 


21.2 


XMM 


2007 Apr 30 


0405380701 


9.8/9.3/2.5 


NGC 7674 


125.3 


XMM 


2004 Jun 2 


0200660101 


8.4/9.2/8.3 



^ Net exposure time after filtering. 
^ For XMM-Newton observations. 

^ Distances based on optical spectroscopic rcdshift using Ho = 70 km s-i Mpc"!, Qm = 0.27 and f7A=0.73 
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Table 2 

APEC model parameters (solar abundance) 



Galaxy 


Nh,i 


kT 


r 


Nh,2 




2pow 


Ipow 




10^2 cm~2 


keV 




1022cm-2 


DOF 


DOF 


DOF 



NGC 0424 


n nt; + 0.04 

0-05-0.03 


0.82+0 j7 


„ Qc+0.32 
/.80_o 28 


16.8l^-j 


269.5 (171) 


273.8 (173) 


846.4 (178) 


NGC 1068 


„ +0.03 

0-31Io.o3 


„ +0.01 

0-61Io.oi 


no + 0.59 

2-02Io.45 


Q 00+1.77 

y-J'J-2.5S 


450.4 (247) 


1013 (249) 


6634 (269) 


NGC 1144 


0.06 


n 07+0.29 
'J-^'-o.oe 


1 r,i+0.37 
1-91-0.24 


47.0^32 


174.7 (149) 


216.8 (151) 


1347 (156) 


NGC 1320 


n n-T+O.OS 


n -70+0. 07 
U.'»_0.07 


"^•"^^-0.19 


'13-5Ii2.3 


269.1 (170) 


311.2 (172) 


639.9 (177) 


NGC 1386^ 


^ ^ .+0.03 

0.04t^^^ 


„ ^^+0.04 


„ „-+0.27 
2-97Io.22 


op. g+19,7 
33-0-13.3 


412.7 (340) 


591.1 (342) 


876.9 (347) 


NGC 1667 


0.05 


00+0.07 


n 1 0+0.34 
^•l*'-0.37 




49.8 (38) 




82.3 (39) 


F05189-2524^'^ 


0.02 


<0.104 


rto+0.13 


^ „c,+0.40 


530.3 (376) 


607.6 (378) 


2212 (379) 


NGC 3982*'^^ 


„ p.„+0.11 
U.03_Q jg 


<0.12 


'J-'^' -0.90 


••• 


21.7 (16) 




45.7 (18) 


NGC 4388 (Chandra) 


-, , +0.49 


<0.18 


„ nn + 0.27 
0-92Io.45 


29.2l4;3 


110.6 (92) 


121.7 (94) 


324.1 (99) 


NGC 4388 (XMM) '° 


0.03 


n on+0.{)l 

0-30+O.01 


-, or: + 0.14 
1-^^-0.10 


26.2+J-2 


580.2 (498) 


844.3 (500) 


3936 (506) 


NGC 4501""^'*'^ 


0.03 


r, .n + 0.16 

0-42Io.o9 


n or\+0.45 

0-30Io.50 




27.5 (38) 


... 


58.9 (40) 


TOLOLO 1238-364^ 


0.06 


^ -o + O.ll 

0-731^.13 


„ ,-+0.31 




51.0 (77) 




72.6 (78) 


NGC 4968^ 


„ „ .+0.10 


<0.13 


-.^+0.41 

1-501^.31 




343.0 (267) 




337.8 (270) 


M-3-34-64 


„ „_+o.oi 


n ^7n+0.02 

0-79Io.o2 


2-68t^.o5 


"4-1 
46.7+1-° 


847.5 (493) 


1660 (495) 


8590 (500) 


NGC 5135^'^ 


0.05 


„ „^+0.24 
0-77Io.22 


^•'^-0.12 


104t»o 


194.8 (132) 


200.8 (134) 


317.8 (138) 


NGC 5194^ 


0.02 


„ ^-+0.05 

0.65+";^^ 


2.20+0 1? 


+62.9 
90.1-43.1 


274.0 (231) 


445.2 (232) 


944.9 (237) 


NGC 5347 


0.02 


<0.24 


-. ,Q+0.24 


63-6+25.7 


31.8 (22) 


36.9 (24) 


78.4 (26) 


iVlrK "iDo 


<U.UD 


r, 70 + 0.03 

0-'3_o 04 


9 09+0.27 


9fi r:+4.9 




DUU.o \^Zvo) 


1 Kr\^ ( o'7f\\ 
louo 


NGC 55068 


11+001 


o.77t«:« 


1 71+0.01 
^•'1-0.01 


2 68+0-03 
^•°°-0.03 


2720 (2385) 


2781 (2387) 


15637 (2389) 


NGC 5506" 


13+001 
'J-i-'-o.oi 


o.85tS:J^ 


1 77+0.01 


2 80+0-01 

^•<5U_o 02 


4171 (3143) 


4299 (3145) 


31991 (3147) 


Arp 220 (XMM)1'5 


0.04 


82+0-05 
^•^^-0.05 


1 97+0.15 




146.0 (145) 




248.4 (147) 


Arp 220 (Chandra) 


47+0-07 














NGC 6890'^ 


<0.10 


0.781°:?^ 


QO+0.8S 
'^•^''-0.74 


+ 18.4 
^'-^-11.3 


164.0 (148) 


171.3 (150) 


197.4 (152) 


IC 5063^0 


64+0-26 


43+0-1^^ 


1 39+0.41 
^•■-"'-0.41 


19.6+2^ 


131.0 (116) 


135.6 (119) 


452.5 (120) 


NGC 7130'^ 


<0.08 


76+0-04 
'^•'°-0.04 


2 41+0-27 
^•^^-0.26 


64.llg-| 


220.7 (199) 


381.8 (201) 


563.9 (206) 


NGC 7172i'6 


0.02 


26+0-03 


1 c:c+0.03 
^•3^-0.01 


7 74+0.09 
'•'^-0.08 


2330 (1748) 


2530 (1750) 


5379 (1751) 


NGC 7582 (XMM)i'5.i2 


0.01 


71+0-01 

^•'^-0.01 


1 95+0-03 
^•=''-'-0.02 


26.011;* 


1586 (886) 


4044 (903) 


15004 (910) 


NGC 7582 (Chandra)^ 


1 24+0.07 
^•^^-0.10 


<0.11 


1 80+0-42 
J^-SU_o 03 


1 Q S+2.3 
I3.0_o 20 


104.9 (81) 


117.5 (83) 


305.2 (85) 


NGC 76741 


0.04 


70+0-13 
^•"^-0.09 


2 92+0-16 
^•^^-0.15 


34.711033 


112.9 (72) 


129.6 (74) 


342.4 (75) 



1 Best-fit Njj was same as Galactic value and therefore frozen at this value. 

2 Best-fit parameters between Chandra and XMM-Newton observations are consistent. 

3 Best-fit parameters between Chandra and XMM-Newton observations are consistent except for the constant multiplicative factor, which 
is much lower for the Chandra observations, indicating extended emission in the XMM field of view. 

4 Best-fit parameters between Chandra and XMM-Newton observations differ due to presence of extended emission in XMM field of view. 
Parameters for the Chandra observation, which isolates the point source, are listed. 

5 Best-fit parameters between Chandra and XMM-Newton observations differ due to variability. 

6 Second power law component normalizations fit indepedently between the two XMM-Newton observations. 

7 Used c-stat. 

* XMM-Newton observations from Feb 2, 2001; Jul 11, 2004; Jul 14, 2004 and Jul 22, 2004. 

9 XMM-Newton observations from Jan 9, 2002; Aug 7, 2004; Jul 27, 2008 and Jan 2, 2009. 

10 Used pileup model. 



28 



Table 3 

Power law model parameters 



Galaxy 


Nh,i 


r 


A^i/,2 






iU cm 




iU cm 




NGC O424I 




2 Q7+0.27 
^•"'-0.26 


16 9+^° 


273.8 (173) 


NGC 4388 (Chandra) 


22+°-2* 


n qo+0.39 
0-^°-0.36 


23.3t3f 


121.7 (94) 


NGC 49682-3 


0.08 


1 Q4+014 




337.8 (270) 


NGC 51352 


0.05 


^•'O_0.10 




200.8 (134) 


NGC 53472 


0.02 


1 41+0-24 


56 2+31-3 


36.9 (24) 


NGC 59532'3'4 


0.03 






39.9 (21) 


NGC 6890^ 


21+°" 




18.9t}6.5 


171.3 (150) 


IC 50632.5 


0.06 


1 48+°-26 

^•^°-0.25 


20.511-1 


135.6 (119) 


NGC 7582 (Chandra)^ 


<0.23 




18.8l2;9 


117.5 (83) 


NGC 76742 


0.04 


2.86lg:}? 


36.91:72^'' 


129.6 (74) 



1 Best-fit parameters between Chandra and XMM-Newton observations are consistent. 

2 Best-fit Nh was same as Galactic value and therefore frozen at this value. 
^ Used c-stat. 

4 Only detected in soft band 

5 Used pileup model. 



Table 4 

2-10 keV X-ray Flux and Luminosity for 12/^m Sample 



Galaxy 


F2-10fceV 


Log L2_l0feeV 


Comments 




10^13 erg/s/cm2 


erg/s 





iNLrO U424 


11. 3 Y 


41.56_(, 


NGC 1068 


54 2+"° 


41 27+0-48 

^l-^'-0.39 


NGC 1144 


00 ^+36.6 
00. i_j^2. 9 


42 77+0-32 
' '-0.21 


NGC 1320 


■H S4+2-83 
■^•^^-1.39 


40 83+°-24 


NGC 1386 


1 55+0-42 


39.48+°- 10 


NGC 1667 


43+0-11 


40.331°;}° 


F05189-2524 




43.00l°-°g 


F08572-I-3915 


<1.26 


<42.02 


NGC 3982 




39.28l°-j2 


NGC 4388 




42 01+°-^* 
'1Z.U1_|3 32 




86.9128-6 


42.08l°;}2 




2441^6 


42 53+°-i2 


NGC 4501 


1 07+°-'^3 


40.17lg:i 


TOLOLO 1238-364 


1 91+0.31 
^•^^-0.26 


40.50+°-}? 


NGC 4968 


2.08l°:2« 


40 65+°-°^ 


M-3-34-64 


32.5+3 1 


42 si+0-04 


NGC 5135 


r, 01 +0.98 

^•■JJ^_1.6g 


40.99l°-jg 


NGC 5194 


1-041^:?^ 


38.95l°-.«5° 


NGC 5347 


2.5811-20 


40.55lg-il 


Mrk 463 


2 95+i-*4 


42 23+°-21 


NGC 5506 




42 75+0-04 
'^-0.05 




11131^5^ 


42 94+0-02 

^^■"^-0.02 


NGC 5953 


<0.51 


<39.73 


Arp 220 


1 n7+°-i8 
i-'J'-o.ie 


40.881°-°? 


NGC 6890 


1 20+4-°i 


40.22o+°5«* 


IC 5063 


i34i;^? 


42 55+0-19 


NGC 7130 


r, r,7+2-09 


41 06+°-^° 

4i.UD_Q 3Q 


NGC 7172 


517l« 


42 96+°-°3 

4Z.9D_o 03 




23411^8 


42 61+0-03 

4Z.D1_Q Q3 


NGC 7582 


21 l+l-'^ 


41 05+°-°^ 




38.6+3? 


41 00+0.03 




164l2f 


41 Qc;+0-42 
4i.90_o 33 


NGC 7590 


<2.72 


<40.17 


NGC 7674 


c- 7-1+3.05 


42.03l°-/J 



Chandra observation 
Chandra observations 

Chandra observation 
Chandra observation 
XMM Jul 2002 observation 
XMM Dec 2002 observation 
Chandra observation 



Chandra observation 

2001 & Jul 2004 observations 

2002, Aug 2004, 2008 & 2009 observations 



2007 observation 
2002 & 2004 observations 
2005 XMM observation 
2001 XMM observation 
Chandra observations 
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Table 5 



Fc Ka Flux and EW 


Galaxy 


Energy 


Global Fit 
cr EW Fluxi 


Energy 


Local Fit 
a EW 


Fluxi 






keV 




keV 





NGC 0424^ 


6, 


45+°-?,:^ 


42+°-" 
^■^^ — 0.11 


4 22+ nn? 

^•^^—0.97 


3 95+°-§? 


6 36+° °! 

^•'-'^—0.05 


21+°-U 
^•^^—0.17 


1 


qq+0.36 

■■''' — 0.39 


9 q7+0.64 
-^.iJ ' —0.69 


NGC 1068^ 


6, 


40+°°? 
^ — 0.01 


<0.03 


65+°-°5, 


c:9+0.41 


6 4+°-°° 

^•^—0.01 


<0.03 


0, 


60+°-°? 

■""-0.05 


5 95+°-iJ 


NGC 1144^ 


6, 


94+0.02 
— 0.02 


<0.07 


0.26+°-°^ 

'-'•■^^—0.06 


1.99+°-*^ 


6.24+°-°? 
^ — 0.02 


<0.07 


0, 


95+0.06 
■ " — 0.06 


1.89+°-if 

^'^— 0.43 


NGC 1320^ 


6, 


07+0.02 
— 0.02 


0.06+°-°^ 

^^^ — 0.03 


3.50+°-i? 


. c:rr+0.22 


6.37+°-°? 

^ ' ' —0.01 


0.05+°-°? 


3, 


.02+°-5?i 
■^^ — 0.50 


1.70+°-?? 

^' ' '-'—0.28 


NGC 1386 (Chandra)^ 


6, 


qq+0.02 
— 0.03 


<0.05 


^ ^— 1.39 


0.51+°?1 

— 0.21 


6.39+°°? 


<0.05 


2, 


on + l.OO 
"^^ — 0.78 


0.72+°-2? 


NGC 16675.6 




6.31 


0.01 




0.16+°?° 

— 0.10 


6.31 


0.01 


0, 


,86+°-?f, 

— 0.50 


0.13+°-?,° 

—0.08 


F05189-25212 (Chandra)^ 




6.14 


0.01 


09+°-°^ 


28+S-?I 

"■^°-0.24 


6.14 


0.01 




<0.17 


<0.59 


NGC 3982 (Chandra)^ 




6.37 


0.01 




0.31+°« 

— 0.22 


6.37 


0.01 






0.11 + ni? 
^ — 0.11 


NGC 4388 (Chandra) 2 


6, 


04+0.02 

— 0.03 


<0.09 


0.31+°-°» 

— 0.08 


4.03+;°6 

— 1.03 


6.34+°°? 

—0.02 


<0.08 


0, 


90+0.11 

— 0.08 


791.47 


NGC 4388 (XMM Jul 2002)^ 


6, 


07+0.01 
— 0.01 


<0.09 


0.46+°-°| 

—0.08 


7.05+;-" 

— 1.16 


6.37+°-°2 

—0.02 


0.08+°-°? 

—0.03 


0, 


.62+°-i° 

— 0.10 


8.69+}-^^ 


NGC 4388 (XMM Dec 2002)^ 






0.06+°-°? 

— 0.02 


0.20+°-°3 

^ ^— 0.03 


9.20+;-2f 

— 1.21 


6.37+°°? 

—0.01 


0.05+°-°? 

—0.02 


0, 


10+0.03 
■ —0.02 


8.39+!-?« 

— 1.14 


NGC 4501 (Chandra)^ 




6.35 


0.01 


<2.28 


<0.29 


6.35 


0.01 




<1.50 


<0.31 


TOLOLO 1238-3643 


6, 


qn+0.29 

— 23 


0.39+°1* 




0-70+030 


6.38+no 


<0.25 


3, 


-17+2.64 
1 96 


0.56+°-^^ 


NGC 4968^ 


6, 


oo+0.'08 

03 




_ 


0.95+2-2? 


f. 07+0.03 

02 


0.07+°-gl 


3, 


0g+o:99 


0.91+°-^^ 


M-3-34-644 


6, 


qn+0.02 

■'-*'-'_0 01 


<0.08 


17+°-04 
^■^ ' _Q 03 


1 9-,+0.32 
' —0 21 


6.30l°:g? 


0.101°:°^ 


0, 


q-t +0.'05 
— 04 


-I 70+0.31 
t. '»-o 23 


NGC 5135^ 


6, 


34+0.04 
■.^^—0.04 


<0.12 


1 18+°-^? 


46+°'?^ 


6 35+°-°t 

U.OU_g 05 


09+°-°! 

"•""—0.06 


2, 


44+0.94 

■ -0.82 


Gl+^il 

^'^^—0.20 


NGC 51942 


6, 


oQ+0.02 
— 0.01 


<0.04 


3.05+°-f, 


0.40+°" 

— 0.08 


6.39+°°? 

—0.01 


<0.05 


4, 


,64+1-12 

— 1.47 


0.57+°-" 

— 0.18 


NGC 53472.6 




6.35 


0.01 


—0.45 


0.37+°-" 

— 0.16 


6.35 


0.01 


1, 


oer+0.53 
— 0.44 


0.45+°-i« 

—0.15 


Mrk 463 (Chandra)^.^ 




6.10 


0.01 


<0.38 


<0.27 


6.10 


0.01 


0, 


-20+°?? 
— 0.13 


0.15+°-i? 

— 0.10 


NGC 5506 (2001 Feb 2)^ 


6, 


00+0.01 

— 0.01 


0.09+°-°? 

— 0.01 


—0.01 


8.02+°-^^ 

— 0.68 


6.38+°-°? 

—0.01 


0.08+°-°? 

—0.02 


0, 


-, 2+0.02 
— 0.02 


7 oc: + 1.15 
— 1.06 


NGC 5506 (2004 Jul 11)2 












6.40+°-°? 

—0.03 


0.20+°-°« 

—0.05 


0, 


-|y + 0.03 

' —0.03 


10.9+?Q 

— 1.9 


NGC 5506 (2004 Jul 14)2 




j: 








6.38+°-°? 

—0.02 


0.10+°-°? 

—0.03 


0, 


-,0+0.02 
' —0.02 


— 1.15 


NGC 5506 (2004 Jul 22)2 






)J 




„ 


6.38+°°? 

—0.02 


o.ii+°-°i 

—0.04 


0, 


14+0.03 

— 0.04 


7.90+1-nK 
— 2.06 


NGC 5506 (2002 Jan 9)2 


6, 


,45+°-°? 

— U.Ol 


0.13+°-°? 

— U.Ol 


o.io+°-°l 

— 0.01 


10.3+°-? 

— 0.7 


6.42+°-°3 

— 0.04 


0.11+°-°? 

— 0.07 


0, 


,09+°-°? 

— 0.03 


yc- + 2.74 


NGC 5506 (2004 Aug 7)2 












6.38+°-°:? 

^^^—0.03 


0.17+°-°! 

" ' —0.04 


0, 


1 9+0.03 
' — 0.02 


11.1+?? 

— 1.9 


NGC 5506 (2008 Jul 27)2 












6.46+°-°2 


0.15t°-°2 


0, 


1 q+0.02 

■ 10 01 


11.7+1-1 


NGC 5506 (2009 Jan 2)2 












6 51+°-°3 
"•'^1-0.02 


n 90+0.05 


0, 


-1 0+0.02 

--'-°-0.02 


is.ol^'l 


Arp 220^ 




6.29 


0.01 


<0.66 


<0.09 


6.29 


0.01 




<0.57 


<0.09 


NGC 6890^.6 




6.35 


0.01 


1 21+l-'"^ 
^•'^^-l.Ol 


'^■^^"-0.13 


6.35 


0.01 


0, 


Qq + 1.28 
-^•'-0.84 


17+°-23 
'^•I'-O.IS 


IC 50633.6 




6.33 


0.01 


05+° °'' 


1 14+0.96 
^■^^-0.94 


6.33 


0.01 


0, 


nt;+0.04 
■'^"-0.04 


58+°-*2 
U.OS_o 42 


NGC 7130^ 


6, 


qn+0.04 

■■''-'-0.04 


<0.09 


70+°-39 


20+°-" 


6 30+° °'' 


<0.10 


0, 


09+0. 48 
■"^-0.33 


26+°-l5 
IJ.^D-O.IO 


NGC 7172 (2007)2 


6, 


qO+0.02 

■''•-'-0.01 


<0.06 


'J-'JO-o.oi 


10+0-56 
■'■^"-0.48 


f. qq+0.02 

"•■'■'-0.02 


ll+°-°'' 
'^•11-0.03 


0, 


in+0.02 
■I'-'-O.Ol 


5 48+1°° 
".^0-0.68 


NGC 7172 (2004)2 


6, 


00+0.01 

■■'"-0.01 


09+° °2 
^■^^-0.02 


12+°-°i 
^•^^-0.01 




f. q7+0.02 
"••"-0.01 


09+° °2 
^•*J^-0.02 


0, 


19+0.01 

■-'-^-o.oi 


q qr\+0.38 
■'•■'"-0.40 


NGC 7172 (2002)2 


6, 


07+0.03 

■■" -0.04 


14+°°6 
^■^^-0.04 


14+°-°3 
^^•^^-0.03 




6 35+°-°4 
"••'^-0.06 


21+°-" 
"•^^-0.06 


0, 


9n+o.06 

■^"-0.04 


^•°1-1.17 


NGC 7582 (XMM 2005)2 


6, 


07+0.01 

■■" -0.01 


<0.04 


41+° °3 
^•^^-0.03 


1 97+0-16 


6.38l°:°i 


n5+°-°i 
^•'JO-o.oi 


0, 


c:o+0.04 
■^^-0.04 


2 40+°-" 
^•*^-0.17 


NGC 7582 (XMM 2001)2 


6, 


07+0.02 

■■" -0.01 


ll+° °3 
'-'■^^-0.02 


62+° °* 
'J-°^-0.07 


q Qq+0.52 

o.»o_o 45 


6.37l°:°i 


<0.07 


0, 


q-, +0.05 

■■'-'--0.05 




NGC 7582 (Chandra)2.6 




6.37 


0.01 


o.i8i°:J°7 


4 07+2.31 
^■■" -1.62 


6.37 


0.01 


0, 


it:+0.07 
■-'-■'-0.07 




NGC 76742.6 




6.22 


0.01 


0.58i°:f7 


49+°-2l 
^■^^-0.23 


6.22 


0.01 


0, 


q7+0.18 
■■" -0.15 


0.48l°:^g 



1 Flux in units of lO—i^crg s~icm~2. Line energies are reported in observed frame. Upper limits on parameters refer to the 90% confidence 
level wrhereas upper limits on the EW and flux signify 3(t error bars. "-" denotes unconstrained parameter. 

2 Fe Ka line detected at greater than the 3a level. 

3 Fe Ka line detected at greater than the 2.5a- level. 
* Fe Ka line detected at greater than the 2a level. 

5 Fe Ka line detected at greater than the 1.5a- level. 

6 XSpec model ZGAUSS used, with E frozen at 6.4 keV (rest-frame) and a frozen at 0.01 keV. 
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Table 6 





Obscuration Diaf 


jnostic Ratios 








Galaxy 


\CJ 1 1 1 \ .obs 


Loj 


-^2— lOfce 


-) 




NGC 0424 


0.22 




0, 


,66 


-2.20 


NGC 1068 


-0.40 




-0, 


.51 


- 


NGC 1144 


1.89 




1, 


,83 


-1.11 


NGC 1320 


0.44 




0, 


,25 


-2.32 


NGC 1386 


-0.71 




-0, 


,72 


-2.63 


NGC 1667 


-0.17 




0, 


18 


-2.81 


F05189-2524 


1.19 




<0, 


53 


-1.53 


F08572+3915 


>2.02 






- 


>-2.52 


NGC 3982 


-0.55 




0, 


14 


-2.63 


NGC 4388 


0.97 




0, 


,45 


-0.94 




1.04 




0, 


,52 


-0.87 




1.48 




0, 


,97 


-0.42 


NGC 4501 


0.46 




0, 


,60 


-2.12 


TOLOLO 1238-364 


-0.58 




-0, 


,09 


-2.99 


NGC 4968 


0.07 




-0, 


,10 


-2.57 


M-3-34-64 


0.33 




0, 


,55 


-1.56 


NGC 5135 


-0.21 




-0, 


,39 


-2.26 


NGC 5194 


0.31 




-0, 


,42 


-2.89 


NGC 5347 


0.77 




0, 


,61 


-2.19 


Mrk 463 


-0.28 




-0, 


,27 


-2.36 


NGC 5506 


1.76 




1, 


,53 


-0.37 




1.94 




1, 


,71 


-0.19 


NGC 5953 


>0.03 




>-0, 


,49 


>-2.91 


Arp 220 


1.77 




<-0, 


,92 


-2.84 


NGC 6890 


-0.29 




0, 


17 


-2.53 


IC 5063 


0.88 




1, 


,14 


-1.03 


NGC 7130 


-0.22 




0, 


12 


-2.33 


NGC 7172 


3.15 




2, 


,17 


0.22 




2.80 




1, 


,83 


-0.13 


NGC 7582 


0.61 




0, 


,38 


-1.38 




0.88 




0, 


,65 


-1.12 




1.50 




1, 


,27 


-0.49 


NGC 7590 


>0.97 




>0, 


,88 


>-1.82 


NGC 7674 


-0.10 




0, 


19 


-2.29 



Table 7 

Correlation of AGN Properties and Star Formation Activity with Obscuration Diagnostics 



Lio/vi vs. 


P 




F2-10fcev/F[O///],o6s 


0.273 


0.053 


F2-10keV /F[oiV] 


0.185 


0.206 


F2-10fceV /Fmih 


0.349 


0.015 


Fe Ka EW 


-0.302 


0.048 


^OIV]/^BH VS. 


F2-10feev/F [0/77], o6s 


0.110 


0.439 


F2-10;se\^/F[o/V'] 


0.056 


0.703 


F2-10keV /FmIR 


0.280 


0.050 


Fe Kq EW 


-0.219 


0.151 


Mbh vs. 


F2-10feeV /FfOJ/I] ,o6s 


0.148 


0.295 


F2-10feev/F[o/V] 


0.115 


0.432 


F2-10fcev/FA//fl 


-0.012 


0.932 


Fe Kq EW 


-0.234 


0.124 


F[o/vi/F[iVei-/i vs. 


F2-10keV /FlOIII],obs 


0.036 


0.804 


F2-10keV /F[oiV] 


0.016 


0.913 


F2-10fceV /FmIR 


0.269 


0.067 


Fe Ka EW 


-0.130 


0.404 


PAW EW 17 ^tm vs. 


F2-10feev/F[o/i-7],o6s 


0.135 


0.346 


F2-10feev/F[O7V'] 


0.062 


0.675 


F2-10keV /F MIR 


-0.055 


0.700 


Fe Ka EW 


-0.192 


0.213 


Q!20-30M'n VS. 


F2-10keV /F[OIII],obs 


0.272 


0.057 


F2-10fcev/F[o/-l,'] 


0.014 


0.922 


F2-10fceV /FmIR 


0.060 


0.675 


Fe Ka EW 


-0.235 


0.127 



Probability that the null hypothesis, 
that the two quantities are uncorrelated, 
is correct. Quantities are statistically 
significantly correlated if P<0.05. 
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APPENDIX 

NOTES ON INDIVIDUAL SOURCES 

NGC 4^4 - The default Chandra aperture extraction spectrum and XMM-Newton spectra were consistent and fit 
simultaneously using a double absorbed power law with a Gaussian component to accommodate the Fe Ka emission; 
including a thermal component did not statistically significantly improve the fit. We note that fitting the 3-8 keV 
continuum with a power law plus a Gaussian led to a tighter constraint on the Fe Ka emission and a more realistic 
EW value. Matt et al. (2003) analyzed the Chandra and XMM-Newton spectra independently using a slightly 
more complicated model, including Gaussian components at 0.55 and 0.90 keV to account for emission features, 
possibly from the OVIII recombination line and the OVIII recombination continuum or Nc IX recombination line, 
respectively. They also added a component for cold reflection, the PEXRAV model in XSpec. However, their 2-10 
keV flux (1.6xl0~^^ erg s~^ cm~^) and Fc Ka EW values (^0.88 keV) are consistent with the values we obtained 
(1.15lg'3Y ^ 10^^^ erg s^^ cm^^ and 1.33lg gg keV, respectively) using a simpler model. 

NGC 1068 - The PN and MOSl XMM-Newton spectra for both observations showed evidence of pileup according 
to epatplot and were therefore not included in the spectral fit; archival Chandra data also exists for NGC 1068, but 
was not used in this study due to the effects of pileup. Also, as several strong emission features were present below 

1 kcV (which are not important for the purposes of this study), we fitted the spectrum from 1 keV to 8 keV. The 
M0S2 spectra were best fit by an absorbed double power law with a thermal model and Gaussian components at 2.0 
keV, 2.43 kcV, 6.4 keV (neutral Fe Ka), 6.66 keV (likely ionized Fc Ka) and 6.95 keV. The neutral Fc Ka line was 
detected at a statistically significant level. Pounds & Vaughan fitted the 3.5 - 15 keV XMM-Newton spectra with two 
continuum components, a cold reflection model (PEXRAV) and a series of Gaussian emission features from 6-8 keV 
(where nine of these features were detected at a significant level). Based on this fit, they find a 3-15 keV flux of 63 
xlO""'^^ erg cm~^, which is consistent with our 2 - 10 keV flux of 54.2l^^°g x lO'^^ erg s ^ cm ^, and a neutral 
Fe Ka EW of 0.60ib0.10 k e V, wh ich agrees with our neutral Fe Ka EW value (0.60±0.05 keV). However, as noted in 
the main text, iMatt et al.l ()2004( ) obtained an EW value of 1.2 kcV, using a PEXRAV model, power law and a series 
of Gaussians for emission line features. 

NGC 1144 - The XMM-Newton spectra were best fit by a double absorbed power law with a thermal component 
and a Gaussian component for the Fe Ka emission. When fitting the 3-8 keV continuum to obtain a local fit 
for the Fe Ka component, a double power law was needed to accommodate the spectrum shape; the power law 
indices of the two components were tied together with the normalizations and an absorption component attenuating 
the second power law allowed to vary. Winter et al. (2008) fit NGC 1144 with the partial covering model in 
XSpec (which is akin to a double absorbed power law model with the photon indices tied together, which we 
have done) and a blackbody component for the soft emission. They derived comparable 2-10 keV flux (3xl0~^^ 
erg s-i cm-2) and Fe Ka EW (0.22 kcV) values as us (33.4t^^| x 10"" erg s'^ cm'^ and 0.25to°t keV, respectively). 

NGC 1320 - The XMM-Newton spectrum to be best fit by an absorbed double power law component with a 
thermal component and a Gaussian component to accommodate the Fe Ka emission. Greenhill et al. (2008) used 
a cold reflection model (PEXRAV) as well as two thermal components (using MEKAL whereas we used APEC) for 
the soft emission; a Gaussian had also been included to model the Fc Ka emission. We derive consistent Fc Ka 
EWs (3.02lg5g kcV vs. 2.20^g'43 keV), yet their 2-10 keV flux is over an order of magnitude higher than ours 
(3.84l^;®^ X 10^1^ erg s^^ cm^^ vs. -43x10"" erg s"^ cm"^). 

NGC 1386 - As noted in the main text, the spectral parameters between the XMM-Newton and Chandra 
observations were consistent, except for a lower multiplicative factor for the Chandra observation, indicating extended 
emission contaminated the XMM-Newton observation. Though the Chandra data were grouped by a minimum of 
5 counts per bin, the XMM-Newton spectra were grouped by a minimum of 15 counts, so we used statistics 
in this analysis. To derive the 2-10 kcV flux, we fitted spectra from both observatories simultaneously to better 
constrain the Chandra spectrum, using a double absorbed power law with a thermal component and a Gaussian 
feature at 6.4 keV to accommodate the Fe Ka emission, yet we report the Chandra flux only. We flt the Chandra 
spectrum independently, both globally and locally, to derive the Fe Ka parameters, binning the data by a minimum 
of three counts and using C-stat. Levenson et al. (2006) fit the Chandra 4-8 keV nuclear spectrum with a 
reflection model (PEXRAV) and a Gaussian component at the Fe Ka energy. We obtain consistent 2-10 keV flux 
values (1.55;°;33 x 10"" erg s'^ cm'^ vs. 2.1±0.1 x 10"" erg s"! cm-^) and Fe Ka EWs (2.30tJ;?° vs. 2.3±1.5 keV). 

NGC 1667 - The XMM-Newton spectra were best flt with a single absorbed power law plus a thermal component. 
To constrain the Fe Ka EW in the local continuum flt (i.e. 3-8 keV), the spectra were binned by a minimum of 

2 counts versus the 15 counts used for the global flt; C-stat was utilized in this local flt. Bianchi et al. (2005) flt 
the spectra with a reflection model (PEXRAV) with a soft excess, including a line at ~0.9 keV. Our Fe Ka EWs 
are consistent (0.86;*Iq'5q vs. <0.60 keV), however our 2-10 keV flux values disagree by about a factor of two 
{0A3toAl X 10"" erg cm"^ g-i vs 10"" erg cm"^ s"^). 
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F05189-2524 - Since the spectral parameters were consistent between the XMM-Newton and Chandra observations, 
other than a lower constant multiplicative factor for the Chandra spectra, we fit these spectra simultaneously to 
constrain the 2-10 keV flux. However, we only report the flux from the Chandra observation as we have demonstrated 
that extended emission contaminates the XMM-Newton field of view. The spectra were best fit by a double absorbed 
power law plus a thermal component, though the temperature was not constrained. The Chandra spectra were fit 
independently to model the Fe Ka emission. Though this feature was not detected, we derived an upper limit on the 
EW of 0.17 keV, consistent with the results of Ptak et al. (2003) who fit the 2002 Chandra observation with a single 
power law plus thermal component (MEKAL). Wc also obtain similar 2-10 keV fluxes {23.5~t^'Q x lO^^""* erg cm~^ 
s~^ vs. 37x10"^^ erg cm^^ s^^). 

NCC 3982 - As noted in the main text, the parameters, fiux and Fe Ka EW values listed are from the model fits to 
Chandra spectrum only as extended emission is present in the XMM-Newton field of view. The spectrum was best-fit 
by an absorbed power law model with a thermal component, using the C-statistic on data grouped by 3 counts. We 
used ZGAUSS to test for the presence of an Fe Ka line, but the EW was unconstrained in both the global and local 
fit, where in the latter, it was necessary to group by 1 count per bin to fit the continuum. Guainazzi ct al. (2005) fit 
the Chandra spectrum in a similar fashion (single absorbed power law with a thermal component) and obtained an 
upper limit on the 2-10 keV X-ray flux of 0.5x10"^'^ erg s""'^ cm~^, which is consistent with our value of O.SG^J'gg 
erg s~^ cm~^. They also report a la detection on the Fe Ka EW of 8±5 keV, which may indicate this parameter is 
unconstrained. 

NCC 4388 - This source varied between the XMM-Newton observations from July to December 2002, increasing 
in flux and decreasing in Fe Ka EW. Both observations were best flt by a double absorbed power law (allowing 
the normalization of the second power law component to vary between the two observations) , a thermal component 
(necessary to flt the soft emission), and a Gaussian component to accommodate the Fc Ka line. The spectral shape of 
the local continuum (3-8 keV, to constrain the Fe Ka EW) for the December observation and Chandra observation 
required a base model of a double power law with an absorption component attenuating the second power law; a 
single power law base model was sufficient to fit the local continuum for the July observation. Beckmann et al. (2004) 
fit the XMM-Newton spectra with a single absorbed power law and a Raymond-Smith thermal plasma model; they 
also detect Fe Ka and a possible Fe K/3 line at ~6.89 keV. Our derived Fe Ka EWs arc consistent (0.62to;io keV vs. 
0.57 keV and 0.18l:^;^^ keV vs. 0.22 keV for the July and December observations, respectively); they do not report 
a 2 - 10 keV fiux or luminosity. The Chandra flux for NGC 4388, from the June 2001 observation, is consistent with 
the XMM-Newton July 2002 flux, though the Fe Ka EW increased, which could be due to the presence of extended 
emission the XMM-Newton fleld of view. Iwasawa et al. (2003) found the nucleus from the Chandra observation to 
be moderately affected by pileup, but we did not see evidence of this when we applied the jdpileup model to the 
spectrum in Sherpa. We obtain consistent Fe Ka EW values as Iwasawa et al. (0.29to io vs. 0.44±0.09 keV), though 
a somewhat higher 2 - 10 keV fiux {74.6^111 x 10"^^ erg s~^ cm'^) than their reported 2 - 7 keV fiux (27x10^^^ erg 
s~^ cm^^), though these values are likely consistent given the error bars on their flux and the more limited energy 
range over which they integrated. 

NCC 4501 - We report the parameters from the Chandra spectral fit as the XMM-Newton observation is 
contaminated by extended emission. The spectrum was best fit by an absorbed power law with a thermal component 
and we utilized the C-statistic as the data were grouped by 3 counts per bin. Brightman & Nandra (2008) also find 
the XMM-Newton field of view to be contaminated by extended emission. They fit the Chandra spectrum with a 
reflection component, ontop of an absorbed thermal power law model. Similar to our work, they do not detect the Fe 
Ka emission line in the global spectral flt. 

TOLOLO 1238-364 - The Chandra spectrum was best fit by an absorbed power law with a thermal component 
and a Gaussian to accommodate t he Fe Ka em i ssion. The data were binned by a minimum of 2 counts and we 
therefore employed the C-statistic. iGhosh et al.l ()2007[ ) fit this spectrum with an absorbed power law plus thermal 
brehmsstrahlung model after binning by a minimum of 20 counts which washes out the Fe Ka feature. They detected 
the line at low signal to noise after re-binning by constant width, but obtain an unconstrained EW whereas we 
detected this feature. 

NCC 4968 - The two XMM-Newton observations for this source were best fit by an absorbed single power law model 
with a Gaussian component at the Fe Ka energy, us ing the C-stat i stic o n data binned by 2 - 3 counts; we saw no 
evidence for variability between the two observations. iBianchi et al.l (j2005| ) fitted these spectra were flt independently 
with a reflection model (PEXRAV) with fixed photon index (F = 1.7), a power law component for the soft excess and 
gaussian for the Fe Ka hue. We obtained consistent 2 - 10 keV flux (2.08j:o 26 xlO^^^ erg cm'^ s"^ vs. 2.7±0.08, 
2.3±0.08 X 10-13 ej.g cm-2 s^i) and Fc Ka EW values (3.06tj^:?^ keV vs. 1.9±0.9, 3.2±1.1 kcV) using the simpler 
power law model. 
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M-3-34-64 - The XMM-Newton spectra were fit by a double absorbed powerlaw with a thermal component and a 
Gaussian at the Fe Ka line energy. Miniutti et al. (2007) fit this source with a reflection model, with the soft emission 
accommodated by a power law model with two thermal components and a photoionizcd gas model and Gaussian 
components at 6.4 and 6.8 keV. Our observed 2-10 keV flux values approximately agree (32.5t3;i x 10~^^ erg s^^ 
cm~^ vs. 21±2 X 10"^"^ erg cm~^), as well as our Fe Ka EWs derived from the global flts (0.171q'o3 kcV vs. 

0. 11±0.02 keV), though our local continuum flt resuhs in a higher EW (0.3llo o4 keV). 

NGC 5135 - Chandra observations of NGC 5135 reveal two X-ray point sources near the nucleus of the galaxy. 
The northern source was identifled by Levenson ct al. (2004) as the active nucleus, so wc restrict our analysis to 
this source, using an extraction region of 1.2." They find the AGN spectrum to be best fit by a model consisting of 
two thermal components, a Gaussian component at ^2 kcV and at the Fe Ka energy, and an absorbed power law. 
We grouped the data by a minimum of 3 counts per bin, employed the C-statistic and find that a double absorbed 
power law with a Gaussian component to accommodate the Fc Ka emission reasonably fits the spectrum. Despite the 
different models used between Levenson et al. (2004) and us, we obtain consistent 2-10 keV fiuxes (2.31ll^° gg x 10~^^ 
erg cm-2 s'^ vs. 2.l0ta[ll ^ 10"^^ erg cm^^ g-i) and Fe Ka EWs (2.44+°;^^ keV vs. 2.4+J j keV). 

NGC 5194 - The nuclear region of NGC 5194 contains several X-ray emitting features: the AGN and diffuse 
emission to the North and South (see Terashima et al. 2001). Chandra is necessary to isolate the Seyfert nucleus and 
we therefore present the results of the Chandra analysis only, and do not include the archival data from XMM-Newton. 
Similar to Terashima et al. (2001), we extracted a source region centered on the optical center of the galaxy with a 
1.5" radius from the Chandra data. The data were binned by a minimum of 3 counts and wc utilized the C-statistic. 
The spectra were best fit by a double absorbed power law with a thermal component and a Gaussian component to 
accommodate the Fe Ka emission. Terashima et al. fit the 2001 observation with an absorbed power law model and a 
reflection model, both of which yield consistent fluxes ('-^1.2x10^^'^ erg s^^ cm" -2) and high EW values {"i.htll keV 
and 4.8^2 5 keV, respectively) which agree with the values we obtain (1.04lQ y3 erg cm' -2 and Amt\fj kcV). 

NCC 5347 - The Chandra spectrum is best fit by a double absorbed power law. To fit the local continuum to test 
for the presence of the Fe Ka line, we rebinned the spectrum by a minimum of 3 counts, utilized the C-statistic and 
detected the line at the 3cr level. Levenson et al. (2006) applied a reflection model (PEXRAV) to the higher energy 
range of the spectrum (4-8 keV) and fit the lower energy portion with power laws, a thermal component and line 
emission. Our 2-10 keV flux and Fc Ka EW values agree, 2.58ti'4i ^ 10~^^ erg s~^ cm"^ vs. 2.2± 0.4 erg s~-^ cm"^ 
and 1.35t!^ keV vs. 1.3±0.5 kcV , respectively. 

Mrk 463 - Extended emission was evident in the XMM-Newton field of view as indicated by the lower observed fiux 
from the Chandra spectrum. Indeed, the XMM-Newton is likely contaminated by the double nucleus (see Bianchi et 
al. 2008), which the Chandra observation is able to resolve. However, other than the constant multiplicative factor, 
the spectral parameters were consistent among the Chandra and three XMM-Newton spectra. The observations were 
consequently fit simultaneously to better constrain the Chandra parameters, though only the flux for the Chandra 
spectrum (for the Eastern source) is reported. To test for the presence of the Fc Ka line, the Chandra spectrum was 
rebinned by a minimum of 3 counts and the C-statistic was employed in the local (3-8 kcV) flt; the line was detected 
at greater than the 99% confidence level according to our simulations. Bianchi et al. (2008) also employed a double 
absorbed power law model to fit the spectrum and we obtain consistent 2-10 keV flux values and Fe Ka EWs, 
2.95^0^2 X 10"" erg s'^ cm'^ vs. 4.1±1.8 xlO'^^ gj.g g-i ^m-^ and 0.20^^:^^ keV vs. 0.2lt[] keV, respectively. 

NGC 5506 - The XMM-Newton MOSl spectra showed evidence of mild pile-up above 6 keV according to the 
SAS tool epatplot for all eight observations and were therefore not flt. For four of these (the 2001, 2002, 2008 and 
2009 observations), the M0S2 spectra was also slightly piled and we excluded them from fitting. Archival Chandra 
data for this source does exist, but were not included in this study because they were severely affected by pile-up. 
We find the spectra to be best fit by a double absorbed power law with a thermal component and several Gaussian 
components to fit Fe K emission, at energies ~6.4 keV (neutral Fe Ka), ^6.7 keV and ~6. 95-7.0 keV; however we 
note that for the local continuum fits, sometimes only two components were needed. The source varies by a factor of 
~1.5 in flux on the time scale of approximately several months, though the Fe Ka EW remains relatively constant. 
For the purposes of our analysis, we use the average flux and Fe Ka EW among the 2001 and July 2004 observations 
(<EW> = 0.14j:°;g^ keV) and among the 2002, August 2004, 2008 and 2009 observations (<EW> = 0.13tg;^^ keV; 

1. e. Figures HI - [TT|) . Guainazzi et al. (2010) studied the Fe Ka emission of these observations in depth, using a suite 
of physically motivated models (i.e. a combination of relativistically/non-relativistically broadened Fe Ka emission 
with relativistic/non-relativistic Compton reflection); the Fe Ka EWs are consistent among these various fits. We 
derive EW values that agree with Guainazzi et al. (2010) using simpler modeling of the local (4-8 keV) spectrum 
with a power law and Gaussian components. 
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Arp 220 - The XMM-Newton and Chandra spectra were best fit by a single absorbed power law model with a 
thermal component; only the absorption varied between the XMM-Newton and Chandra observations, though this 
had a negligible impact on the observed flux. We discarded the XMM-Newton spectra from 2005 from our fitting due 
to low signal-to-noise, though we note the best-fit parameters were consistent with the other XMM-Newton spectra. 
To test for the presence of Fe Ka emission in the nuclear region, we utilized the Chandra data only and rebinned by 
a minimum of 3 counts, employing the C-statistic, but the line was no t detected. We note, however, that an emission 
line for ionized Fe Ka at E = 6.51 keV was detected, consistent with llwasawa et all ()2005[ ). The Chandra data were 
first analyzed by Clements et al. (2002) who report a double nucleus and a halo of extended emission. We obtain 
consistent fiuxes between their 3" extraction area (which encompasses the double nuclei) and our 4.5" extraction 
region: 1.07^° ie ^ 10~^^ erg s~^ cm~^ vs 1.0x10"^^ erg s~^ cm~^ 

NGC 6890 - The XMM-Newton spectra were grouped by a minimum of 3 counts per bin and were fit wit the 
C-statistic. The data were best fit by a double absorbed power law model and the Fe Ka line was marginally detected 
at the '^93% confidence level. Shu et al. (2008) fit this source with a single power law and did not detect the Fe 
Ka line, though this likely due to their choice of binning the data by a minimum of 20 counts which would eradicate 
the weak Fe Ka feature. We obtain consistent 2-10 keV flux values given the error range on our derived flux: 
l.2Qtt^l\ ^ 10"" *3i'g cm-2 vs. 0.69x10^13 erg s^^ cm^^. 

IC 5063 - The Chandra spectrum was moderately affected by pileup, ^14% according to the jdpilcup model in 
Sherpa. We therefore included a pileup model in the XSpec spectral fits, allowing only the grade migration parameter 
(a) to be free, and obtained an a value of 0.37; excluded the pileup component when calculating the observed 2-10 
keV flux. The jdpileup model indicated no pileup in the local 4-8 keV spectral fit, so it was modeled without a pileup 
component in XSpec. The broadband spectrum was best fit by a double absorbed power law and our simulations 
indicate that the Fe Ka emission feature is significant at greater than the 2.5a- level. This Chandra spectrum has not 
been previously analyzed. 

NGC 7130 - The Chandra spectrum of NGC 7130 was best fit by a double absorbed power law with a ther- 
mal component for the soft emission and a Gaussian feature at the Fe Ka energy; we used the C-statistic with 
the data binned by 2 counts. This source was studied in detail by Levenson et al. (2005) where they fit the 
AGN spectrum with a double absorbed power law with a gaussian component as well as two thermal compo- 
nents. Though they added an extra thermal component, our derived 2-10 keV flux and Fe Ka EW values are 
consistent {2.07tlM x 10"" erg cm'^ s'^ vs. l.Gtol x 10~" erg cm'^ s'^ and 0.&2toil keV vs 1.8^!^:^, respectively). 

NGC 7172 - The XMM-Newton observations were best fit by a double absorbed power law with a thermal 
component and Gaussian feature at the Fe Ka energy. The normalization of the second power law was consistent 
between the the 2002 and 2004 observations, yet had to be fit independently for the 2007 observation. We only fit 
the FN and M0S2 spectra for the 2007 observation as the MOSl spectrum showed evidence for milder pileup at 
higher energies from the task epatplot. Our results indicate that the source increased by about a factor of 2 in flux 
between 2004 and 2007. The Fe Ka emission features were fit independently among the three observations, though 
we use the average Fe Ka EW for the 2002 and 2004 observations in the plots (<EW>=0.16;^;!^^ keV; i.e. Figures 
m - [Tl]) as the 2-10 keV flux values are consistent between the two observations and the EW values are not widely 
discrepant in both the local and global fits, indicating that the variations in the independent Gaussian fits are likely 
not significant. Noguchi et al. (2009) fit the 2007 observation with a double absorbed power law, thermal component 
(using MEKAL whereas we used APEC) and two Gaussian components, one at the Fe Ka energy and the other at 
1.7 keV; we obtain consistent 2-10 keV fluxes (517^40 x 10"^^ erg cm~^ vs. 423x10"-^^ erg cm~^ s"-^) and Fe 

Ka EWs (0.10j:!^:oi keV vs. 0.07±0.01 keV). Shu et al. fit the spectra from the 2002 XMM-Newton observation by 
a double absorbed power law with a Gaussian component at the Fe Ka energy, similar to our analysis, though they 
do not include a thermal component. Our 2-10 keV flux values are consistent (234^i[g x 10"^^ erg cm~^ s"^ vs. 
220x10"" erg cm'^ s"^), but our Fe Ka EWs are not (0.14;o o3 keV from the global continuum fit vs. 0.04±0.03 
keV). This discrepancy could be due to constraints placed on their modeling of the Fe Ka line: they froze the energy 
at 6.4 keV, whereas we allowed this parameter to be free and it is not clear whether they placed a similar con- 
straint on (T. Regardless of this discrepancy, both EW values are consistent with a Compton-thin source. Analysis of 
the EPIC data for the 2004 observation has not been previously published, making our analysis the first for this datasct. 

NGC 7582 - This source varied between the 2000 Chandra observations and later XMM-Newton observations, as 
well as between the 2001 and 2005 XMM-Newton observations; the Chandra and XMM-Newton observations were 
fit independently. The Chandra spectra were moderately affected by pileup (^^30 - 49% according to the jdpileup 
model in Sherpa) and were therefore fit with a pileup model component in XSpec, with only a, the grade migration 
parameter, allowed to vary; the pileup component was discarded before calculating the fiux and the Fe Ka EW. 
The broad-band Chandra spectra were best fit by a double absorbed power law whereas the XMM-Newton spectra 
required a thermal component and Gaussian components at the Fe Ka energy and at ionized Fe Ka energies (^^6.72 
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keV for the 2005 observation and ~6.97 keV for the 2001 observation^], ~1.84 kcV and ~2.47 kcV. The Fe Ka feature 
was detected at a statistically significant level for all observations according to our simulations. NGC 7582 dimmed 
between the Chandra observation and each subsequent XMM-Newton observation and the Fe Ka EW increased. This 
decrease in flux with increase in Fe Ka EW could reflect a variation in the obscuring medium, where the obscuration 
enhanced over time. Indeed, such an interpretation is favored by Piconelli et al. (2007), who postulate the existence 
of multiple absorption components in this system: a higher column-density absorber (possibly mildly Compton-thick) 
attributed to the putative torus and a lower-column density absorber acting as a screen to both the reflected and 
transmitted radiation. Though Piconelli et al. (2007) fit the XMM-Newton observations with a more complex model 
(PEXRAV) we obtain consistent fluxes (21.lti;g x IQ-^^ erg s^^ cm"^ vs. 23.5x10-13 erg s^^ cm'^ for the 2005 
observation and 38.6l3 ° x lO^^"^ erg s^"'^ cm^^ vs. 40.2 xlO^^'^ erg s^^ cm^^ for the 2001 observation), though lower 

Fe Ka EW values (2005: 0.58l°;°4 keV vs. 0.77^°;°^ keV; 2001: 0.31^°;°^ keV vs 0.62^°;°^). Dong et al. (2004) lit 
the Chandra spectra independently with a double absorbed power law, yet obtain an observed flux about half of ours 
(164^87'^ X IQ-i'^ erg s^^ cm^^ vs. ^75 xlO^^'^ erg s'^ cm^^), though this discrepancy could result from us modeling 
pilcup whereas they did not; we obtain consistent Fe Ka EW values (0.151q'o7 keV vs. an averaged ~0.19 keV). 

NCC 7674 - The global XMM-Newton spectra were best fit with a double absorbed power law. To fit the local 
continuum between 3-8 keV, the data were binned to a minimum of 3 counts and the C-statistic was utilized; the 
Fe Ka feature was detected a statistically significant level. Analysis of the broad-band XMM-Newton spectra for this 
source has not been previously published. 

SIMULATED NUSTAR DETECTIONS: 10-40 KEV 

Using ARE, RMF and background files provided by the NuSTAR team, which include separate ARE and background 
files for a 45" and 101" point spread function (useful for weak and strong sources, respectively) we generated a simulated 
NuSTAR spectrum in the 10-40 keV energy range as described in the main text. If the net count rate was > 10"^ 
s""'^, we utilized the simulated spectrum corresponding to the larger PSE. otherwise we used the spectrum generated 
with the smaller PSE. In Tables [BT] and lB2l we list the simulated source count rate and corresponding exposure time 
for the target to be detected at greater than the 5a level above the background, which was either 8 x 10-"*^ s-^ or 
^ 4 X IQ-"^ s-i, depending on the PSE. For the sources where this derived exposure time is under 5 ks, we instead 
list the exposure time necessary for at least 100 counts to be detected; if this exposure time is also under 5 ks, we 
adopt a minimum exposure time of 5 ks. For the three 2-10 keV non-detections from the 12/im sample, we list the 
exposure times using a spectrum simulated from the PLCABS model in XSpec, with Nh = 10^^ cm^^, the number 
of scatterings set to 5, r=1.8 and the normalization adjusted such that the 2-10 keV flux equals the 3ct upper limit. 
Again, we note that such a flux estimate is quite optimistic and the corresponding derived exposure times necessary 
for detection should be considered lower limits and are listed as such in Table IbTI 



7 However, in the local (3 - 8 keV) fits, this higher energy Gaus- ^ ^ ^j^^ 2^^^ observation, and shifts to a best-fit 

sian IS more consistent with te Kp emission, with a best-fit ccntroid i j r c m i r j-i, onoc: u i- 

^ ' centroid energy oi b.yl keV lor the 2005 observation. 
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Table Bl 

NuSTAR Simulation Summary for 12^m Sample 



Galaxy Simulated Source Count Rate Exposure Time^ 





counts/sec 


ks 


NGC 0424 


3.48x10"^ 


8.8 


NGC 1068^ 


6.53x10"^ 


5.0 


NGC 1144^ 


1.77x10"^ 


5.0 


NGC 1320 


1.33x10"^ 


30 


F05189-2524^ 


3.58x10"^ 


5.0 


F08572+3915^ 


1.10x10"^ 


>9.1 


NGC 3982 


1.19x10"^ 


35 


NGC 4388 


9.60x10"^ 


5.0 


NGC 4501 


5.04x10"^ 


5.7 


M-3-34-64^ 


1.08x10"^ 


5.0 


NGC 5135^ 


6.88x10"^ 


15 


NGC 5194 


1.82x10^^ 


20 


NGC 5347^ 


1.77x10"^ 


5.6 


Mrk 463 


3.90x10"^ 


7.7 


NGC 5506^ 


1.51 


5.0 


NGC 5953 


3.72x10"^ 


>8.2 


Arp 220 


9.00x10""* 


53 


IC 5063^ 


6.22x10"* 


5.0 


NGC 7130 


3.14x10"^ 


10 


NGC 7172^ 


7.67x10"* 


5.0 


NGC 7582^ 


4.68x10"^ 


5.0 


NGC 7590^ 


3.52x10"^ 


>5.0 


NGC 7674^ 


6.17x10"^ 


16 



Exposure time needed for source to be detected above the background at 
greater than the 5(t level. 

Exposure time listed is for a detection of 100 counts as a 5(t detection 
above the background is <5 ks. 

Minimum exposure time of 5.0 ks is adopted as the exposure times for both a 5a" 
detection above the background and for a detection of at least 100 counts are <5 ks. 



Table B2 

NuSTAR Simulation Summary for the [OIII] Sample 
Galaxy Simulated Source Count Rate Exposure Time* 





counts/sec 


ks 


NGC 0291 


1.40x10"^ 


28 


Mrk 0609^ 


1.72x10"^ 


5.8 


IC 0486^ 


1.27x10"* 


5.0 


2MASX J08244333+2959238^ 


7.36x10"^ 


14 


CGCG 064-017^ 


1.12x10"^ 


8.9 


2MASX J11110693+0228477 


1.68x10"^ 


22 


CGCG 242-028^ 


1.22x10"^ 


8.2 


SBS 1133+572 


1.14x10"^ 


37 


Mrk 1457^ 


6.84x10"^ 


15 


2MASX J11570483+5249036^ 


5.87x10"=* 


17 


2MASX J12183945+4706275 


3.23x10"^ 


9.7 


2MASX J12384342+0927362^ 


2.31x10"^ 


5.0 


CGCG 218-007^ 


1.34x10-2 


7.5 



Exposure time needed for source to be detected above the background at greater than 
the 5(T level. 

2 Exposure time listed is for a detection of 100 counts as a 5cr detection above 
the background is <5 ks. 

^Minimum exposure time of 5.0 ks is adopted as the exposure times for both a Scr 
detection above the background and for a detection of at least 100 counts are <5 ks. 



